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precipitation imbedded in the summer monsoon trough over 
East Asla. This work used a two dimensional’ frontogenesis 


model to study the Mei-Yu front, especially its associated 
low-level jet (LLJ) which ee oe cge during inten 
convection. Several numerical experiments were carried out 
to simulate the quasi-steady frontal structures in different 
environments. The results resemble many observed features, 
such as the stronger temperature and moisture gradients in 
the midlatitude fronts and the stronger horizontal shear and 
LLJ in the subtropical fronts. The LLJ is most conspicuous 
in the subtropical simulation when high humidity and surface 
fluxes are included. The results further suggest that the 
LLJ is developed through the Coriolis torque that is exerted 


by the low-level poleward branch of a "reversed Hadley cell" 


equatorward of the front. This thermally direct cell is 
different from the normal cross~frontal secondary 
Cie eu lata one Its development depends on the occurrence of 
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1a INTRODUCTION 


A review of recent observational studies of the main 
phenomena of this study, the Mei-Yu, and the associated 
disturbances will be presented first. There follows a 


review of the pertinent theoretical studies. A description 


of the objectives of this numerical study of the Mez=iiae 


front will then be given. 


A. OBSERVATIONAL STUDIES 
Be et 

Over eastern ASia, the transition period from the 
northeast winter monsoon to the southwest summer nonseam 
possesses a distinctive climatological phenomenon which is 
important for regional weather forecasting. The seadaul 
rainfall reaches a maximum during this period (Fig. 1.1). 
The rainfall during this period is called "Mei-Yu" (plum 
rain) in China and "Baiu"” in Japan. During this period, the 
rainfall may be continuous or intermittent for several days 


to a few weeks and is characterized by frequent rain-showers 


and thunderstorms, with rainfall rates of up to ae few 
hundred millimeters per day (e.g., Chen, 1983; Ninomiya, 
1980). The Mei-Yu season first appears in southern China in 


May, and influences the weather over southern China and 


Taiwan from mid-May to mid-June (Chen and Chi, 1980). ie 


eZ 


then migrates over the Yangtze Valley and Japan from mid- 
Hames to omad=J u lye 

Synoptically, the Mei-Yu rainfall is associated with 
the repeated occurrence of a front that develops in the 
midlatitudes and slowly moves southeastward to a quasi-~ 


stationary position extending from southern Japan to 


southern China. Satellite pictures” reveal a nearly. 
continuous cloud band accompanying this front, which is 
often called a Mei~Yu front or a Baiu front. The frontal 


position coincides with the early summer monsoon trough in 
eastern Asia and its existence may be viewed = as a 
manifestation of the trough. The beginning of the rainy 
season also signifies the onset of the northern summer 
Monsoon in this part of the world. The Mei-Yu often 
undergoes sequential jumps between May and July (Fig. 1.2; 
Tao and Chen, 1987) depending upon the strength and the 
position of the North Pacific Ocean subtropical high (Tao 
endeDingss 1981). As the Mei-Yu onset occurs in the Yangtze 
River region in June, the Indian monsoon a also being 
established and the westerly jet stream begins to retreat to 
the north of the Tibetan Plateau. At the same time, the 
subtropical high in the western Pacific area moves northward 
poe cO0>=25en The Mei-Yu season lasts through mid-July, when 
the subtropical high moves northward to 30°N. At that time, 
the rain belt is to the yams of the Yellow River, which 


Signifies the beginning of heavy rainfall period in northern 


China and the ending of the Mei-Yu in subtropical China (Tao 
and Chen, 1987). 
2. Mei-Yu Front 
The climatology of the Mei-Yu front over southern 
China and the adjacent areas was examined by Chen and Chi 


(1980) for the periods of 15 May-15 June, 1968-1977. They 


found that the Mei-Yu frontogenesis takes place in the 


subtropical part of southern China to the south of ~35°N. 
All the Mei-Yu fronts were formed in the area of 20°-35°9N 
and 100°-l130°E (Fig. 1.3). The lifetime of individual cases 
of Mei-Yu fronts ranges from 3 days to 22 days with an 
average of 8 days. The kinematic background for the 
frontogenesis 1s the deformation field between the 
subtropical high over the Pacific amd the migratory Sigg 
over central China. They suggested that the sea-surface 
temperature contrast in the nearby South and East China Seas 
may also be a favorable factor for the frontogenesis. 

In contrast to the polar front, the temperature 
contrast across the Mei-Yu front is typically small, but 
there is always a marked moisture gradient near the frontal 
zone (Akiyama, 1973a; Chen and Tsay, 1978; Matsumoto et-al., 
1970, 1971; Tac,  IS80)- An analysis of the kinematic 
structure of the Mei-Yu frontal system over Southeast Asia 
was carried out by Chen (1978). The Mei-Yu frontal system 
possesses a rear aoe shear line which roughly coincides 


with a trough at 850 and 700 mb, although the associated 
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geopotential gradient is rather weak in the immediate 
Prcinity wa trowgen. The shear line approximately coincides 
with the maximum north-south gradients of mixing ratio and 
equivalent potential temperature. To the south of the shear 
iene, ariowelevelmjcte(LidJjsthat ws characterized by maximum 
mixing ratio and high equivalent potential temperature is 
aad near ’0Ommbeet rae. 2.4).  Sitrenge cyclonic vorticity, 
horizontal convergence and upward motion are observed along 
ead ae = the south of ,the shear line. The cloud band 
approximately coincides with the area of maximum cyclonic 
Merticity , horizontal convergence and upward motion. 
Intense convective clouds in the southern portion of the 
cloud band tend to be associated with the maximum values of 
these kinematic parameters as well with relative humidities 
greater than 80% in the lower-middle troposphere. Clear 
Miares Ower the East Chima’ Sea sand mts vicinity are 
accompanied by downward motion of dry air (relative humidity 
less than 30% at 850 mb, less than 10% at 700 and 500 mb). 
Chen and Chang (1980) studied the structure = and 
calculated a vorticity budget for a Mei-Yu system during the 
meriod WO-%5 Jue 1975. During the mature and decaying 
stages of the monsoon trough, the structure of the eastern 
section (southern Japan) and the central section (East China 
Sea) is different from that of the western section 
(southeastern China). In the eastern and central section, 


the structure resembles a typical midlatitude baroclinic 


front with a strong vertical tilt toward an upper-level cold 
core low and a strong horizontal et. - veel gradient. On 
the other hand, the western section resembles a tropical 
disturbance with an equivalent barotimepae, warm core 
structure, a relatively weak horizontal temperature gradient 
and a rather strong horizgontal wind shear in the lower 
troposphere. In the mature stage, the upward motion regal 
maximum values at 500 mb. It is located over the &50 mb 
monsoon trough in the western (W) section, 4° latitude south 
of the trough in the central (C) and 3° latitude southeast 
of the trough in the eastern (E) section. On the north side 
of the trough, downward motion prevails with a maximum at 
400 mb in sections C and E, and at 700 mb in section W. 
This well-organized vertical velocity couplet is: indicative 
of a well-developed secondary circulation normal to the 
trough. In the decaying stage, the circulation becomes 
substantially less vigorous in sections W and C, ‘but it 
maintains intensity in section E. Tsay and Chen (1980) 
analyzed the processes which force the vertical motion in 
sections W, C and E. They found the eastern portion 
partially maintains the characteristics of a midlatitude 
system, while the western portion behaves similar to a 
tropical system. 

The Mei-Yu front can be several thousand kilometers 
long, and therefore et is a large-scale monsoonal 


phenomenon. However, the cross~front scale is considerably 
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smaller, which is within the mesoscale range. Among the 
many mesoscale features that are observed .near the Mei-Yu 
frontal zone, the most important are the organized mesoscale 
convective systems embedded within the broad Mei-Yu cloud 
mend whaech liesWon the warm side of the front (Fig. 1.5). 
These convective systems tend to move along the front from 
west to east and produce locally heavy rainfall rates of | 
> 100 mm day-1 (Akiyama, 1974, 1978; Chen, 1977; Ninomiya 
eee Akiyama, 19%, W972, 1973, 1974). 
3. Low Level Jet (LLJ) 

To the south of the Mei-Yu (or Baiu) front, there is 
Omeen aestromg ageostrophic low level jet Chiiiimewhichiaies 
usually observed between 1.5 and 3 km (850-700 mb) and is 
strongly correlated to the heavy rainfall (Akiyama, 197%3b; 
Chen, 1977, 1986; Matsumoto, 1972; Matsumoto, et al., 1971; 
Ninomiya and Akiyama, 1974; Tao and Chen, 1987). Akiyama 
(1973a) observed that the heavy rainfalls in the Japan Baiu 
season occur most frequently to the north of the low level 
jet stream and seldom occur in areas more than 200 km south 
of the jet axis. Matsumoto (1973) found that the peak of 
heavy rainfall precedes by 6 h the peak of the low level jet 
speed. He also found that the peak of heavy rainfall is in 
phase with the maximum acceleration of the LLJ. 

In the Great Plains and Midwest areas of the United 
States, the low level jets are a climatological feature with 


the maximum frequency during the spring and summer months 


iL 


(Bonner, l1968)> The low level jets “in hs area Sane 


basically a boundary layer phenomenon that are characterized 


by a diurnal oscillation. The jet reaches maximum intensity 
by early morning, and is associated with a nocturnal 
temperature inversion (e.g., Blackadar, 1957; Bonner, 1968; 
Holton, 1967; Lettau, 1967; Wexler, 1961). Different 


theories (discussed later in section B) have been proposed. 
to explain the generation of the low level jet and the 
seasonal, temporal and geographic variations. In the British 
Isles, a low level jet in conjunction with midlatitude 
cyclones was studied by Browning and Harrold (1970) using 
Doppler radar observations. Browning and Pardoe (1973) 
found that for six cases of ana-cold fronts occurring in the 
colder months in British Isles, the low level jet had a 
maximum velocity in the layer 900 to 850 mb of 25 to 30 m gs} 
and occurred during night or day. 
4. Mesoscale convective systems (MCS’s) 

Rainfall characteristics over Taiwan and vicinity 
were studied by Chen (1977) and Chen and Tsay (1978) using 
dense hourly rainfall reports, radar echoes and satellite 
pictures. The low level southwesterlies, especially the 
LLJ, tends to create potential instability to the warm side 
of Mei-Yu front. Continued large-scale ascent then leads to 
the release of potential instability through organized 
mesoscale convective systems. The mesoscale convective 


systems embedded within the broad Mei-~Yu cloud band can be 
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or several different sizes (meso -a, -8, or -—y). A case 
study by Chen and Tsay (1978) revealed the existence of 
these mesosystems in Taiwan with wavelengths of 60, 85 and 
DaGekmneand the corresponding periods of 2.6, 3.5 and 7.1 h. 
Peecpectrum analysis @Taoy 1980) of radar echoes in the 
rainbands along the Mei-Yu front over the Chinese mainland 
also showed the existence of organized mesoscale convective 
systems with wavelengths of 150, 84 and 50 km. Akiyama 
(1978) found that the echoes within the rainband in the Baiu 
season in Japan are organized into mesoscale clusters with 
spatial intervals on the order of 100-150 km. 

Simulations of these mesoscale convective systems 
require a three-dimensional model. In our two-dimensional 


Simulations, we will not study these phenomena. 


Mm.  CREORETICAL STUDIES 
i. eaten. t 

Quasi-geostrophic frontal theories (Stone, 1966; 
Williams, 1968; Williams and Plotkin, 1968) showed that 
large-scale deformation fields can lead to surface 
frontogenesis and result in strong horizontal temperature 
gradients that are characteristic of fronts. However, the 
structure of the front and the associated secondary 
circulation do not tilt with height. It also lacks cyclonic 
wind shear at the front and has several other unrealistic 


features, including the infinite time required for the front 


ie 


to form. Later studies using primitive equations (Williams, 
1972) and the semi-geostrophic equations (Hoskins and 
Bretherton, 1972) removed these unrealistic features by 
including the ageostrophic advections of temperature and 
momentum. In particular, Hoskins and Bretherton (1972) used 
the set of semi-geostrophic equations to show that the 
nonlinear effects of the ageostrophic advections can be 
represented by the quasi-geostrophic solutions, if these 
solutions are written in the geostrophic coordinates and 
then transformed to the physical space. Later, the semi- 
geostrophic theory was expanded to three dimensions by 
Hoskins (1975) and Hoskinsvand Digaanwe. (197752 

Frictional and diabatic effects also have _ been 
studied by a number of investigators. Williams (1974) found 
that vertical mixing can determine the slope of a front, 
and, together with horizontal mixing, can determine the 
cross-width of the frontal zone. Keyser and Anthes (1982) 
showed that the proper simulation of the concentrated upward 
motion in the surface frontal zone depends on the planetary 
boundary layer parameterization. Williams et al. (1981), 
Ross and Orlanski (1978), and Hsie and Anthes (1984) found 
that latent heat release intensifies both the upper level 
westerlies and the low level easterlies poleward of the 
surface front, and the ageostrophic secondary circulation 


across the front. 





2. Symmetric Instability 

Within extratropical cyclones, banded structures in 
weiid's and precipitation have been observed boetn in 
conditionally stable and unstable Wes Teed: (Hobbs, 1978). 
Bennetts and Hoskins (1979) suggested that these bands are a 
result of conditional symmetric instability. Emanuel (1979) 
demonstrated that Clureul ations resulting froin tnaS 
instability are fundamentally mesoscale at character. 
Recent observations indicate that conditional symmetric 
instability is responsible for some precipitation bands 
within extratropical cyclones (Bennetts and Sharp, 1982; 
Moore, 1985; Parsons and Hobbs, 1983; Sanders and Bosart, 
FSeo") . 

Bennetts and Hoskins (1979) -and Emanuel (1982) 
assessed the degree of conditional symmetric instability 
based on infinitesimal displacements in saturated 
atmospheres. Emanuel (1983a,b) considered parcel 
instability which leads to a measure of the stability of the 
moist baroclinic atmosphere to finite slantwise reversible 
displacements of a two-dimensional air parcel or air "tube". 
Dynamically, inertial, symmetric and convective 
instabilities are very closely related. Each of these 
instabilities may be thought of as resulting from an 
imbalance of body forces acting on a fluid parcel. In the 
case of convective instability, the gravitational force acts 


vertically. In the case of inertial instability, the 


om 


centrifugal force acts horizontally. Symmetric instability 
simply represents motion resulting from a combination of 
these forces, and thus the motion is "slantwise.”" Following 
Emanuel (1983a,b), iat _ be shown that slantwise 
instability requires the slopes of the pseudo angular 
momentum (M) surfaces to be smaller than those of the 
equivalent potential temperature (Q90e) surfaces. This. is” 
equivalent to stating that @e decreases along the constant 
M surfaces upward or that M decreases in the cross-front 
direction toward warm air along a 6e-~surface. From the 
definition of M, the latter implies that the absolute 
vorticity is negative along 9 e-surfaces. 

A case study by Emanuel (1983b) of ae strongly 
baroclinic .flow with observed linear bands of clouds and 
precipitation reveals that the atmosphere is initially 
conditionally unstable to upright convection and tends to 
become neutral with respect to slantwise displacements. 
Such an atmosphere possesses moist adiabatic lapse rates 
along M surfaces, even though it may be substantially stable 
to vertical displacements. In this case, conditionally 
unstable baroclinic flow undergoes a form of slantwise moist 
adiabatic adjustment that renders the atmosphere neutral to 
slantwise reversible displacements. The adjustment uses up 
the total slantwise convective available potential energy 
(SCAPE) rather than just the convective (gravitational part 


of ) available potential energy (CAPE) for upright 





Convection. After such an adjustment, the SCAPE of the 
atmosphere is very close to zero while the CAPE actually 
becomes negative. 

Emanuel (1985) carried out a diagnostic study using a 
semi-geostrophic model with small potential vorticity (1l.e., 
small symmetric stability) in the region of upward motion 
mae not so small potential vorticity in the region of. 
downward motion. He demonstrated that frontal forcing 
mesults in a strong, concentrated, sloping updraft that is 
located § ahead of the region of maximum gepstrophic 
frontogenetical forcing. The distance between the "front" 
and the updraft will, in the atmosphere, have a typical 
value on the order of 50-200 km. This was observed by. 
Sanders and_ Bosart CtSeSai, who found that slantwise 
convection can persist even after a condition of near- 
neutrality has been achieved under frontogenetical forcing. 
It has also been shown by Thorpe and Emanuel (1985) that in 
an atmosphere which has small stability to slantwise 
convection, frontogenesis proceeds at a greatly increased 
rate (due to latent heat release) and the horizontal scale 
of the ascent in the warm air is considerably reduced. 
eross-iromtals circulations produced by geostrophic forcing 
are similar in structure to slantwise convection, despite 


the (small) stability to the latter. 


2s 


3. Low Level Jet 2m) 

In the Vie Wiey, of Japan, the Balu frome 1s 
accompanied by an ageostrophic LLJ near 700 mb, which is 
oriented along the heavy rainfall zone. The ageostrophic 
LLJ is suggested to be a result of the downward transport of 
horizontal momentum caused by cumulus convection (Akiyama, 
LUT 3sa. Matsumoto, T9372: Matsumoto and Ninomiya, ] Sige 
Matsumoto et al., 1971; Ninomiya, 1971; Ninomiya and 
Akiyama, 1974). In China, several different relationships 
between the low level jet and heavy rains have been observed 
in the Mei-Yu season. In one of these, a LLJ is coupled and 
parallel to the upper tropospheric jet. The formation and 
strenthening of the LLJ is often associated with the heavy 
Pownall 1... Chen (1982) suggests that coupling of upper and 
lower level jets can be explained by an unstable inertia- 
gravity wave that arises from a thermal wind adjustment 
process that corrects an imbalance due to the advection in 
the entrance region of the upper level jet. This coupling 
is in contrast to that studied by Uccellini and Johnson 
(1979) in the Great Plains of the United States, where the 
LLJ tends to be perpendicular to and in the exit region of 
the upper level jet. 

several theories have also been proposed to explain 
the formation and diurnal oscillation of the LLJ ine 
Great Plains region of the United States. Blackadar (1957) 


suggested that oscillation in the boundary-layer mixing 
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processes is responsible ™for the diurnal variation of the 
LLJ intensity which reaches maximum values by the early 
morning. The temporal and geographic preference of the LLJ 
— been further explained by diurnal radiation fluxes over 
sloped terrain (Holton, 1967; Hsueh, 1970; Lettau, 1967), 
and topographical characteristics (Paegle and Rasch, 1973; 
Paegle, 1978). Observations also showed that low level jets 
may also develop in response to synoptic or subsynoptic 
processes, particularly through a response to leeside 
cyclogenesis that is common to the Great Plains (Djuric and 
Damiani, 1980; Djuric and Ludwig, 1983; Naistat and Young, 
eres Newton, 1956; 19867; Reiter, 1969; Uccellini, 1980). 
Numerical simulations and a case study by Uccellini and 
Johnson (1979) suggested that Rema Gucilepiment of the LLJ is 
coupled to the upper tropospheric jet streak by a mass 


adjustment process within the exit region of the streak. 


Cm OontECTIVE OF THIS STUDY 

The objective of this study is to use a two-dimensional . 
frontogenesis model similar to that of Williams et al. 
(1981) to study the quasi-steady state structure of the Mei- 
Mietront sand the LLJ. The model contains a ie ve~ 
adjustment type cumulus heating parameterization and a bulk 
aerodynamics formula for the boundary layer. To study the 
possibility of cumulus momentum transport or the diurnal 


oscillation processes, additional parameterizations would be 


necessary. Since considerable uncertaintwes exist tinesbomh 
the cumulus momentum transport and the diurnal radiation 
parameterizations, it is desirable to carry out a modeling 
study excluding these processes to see if the LLJ may be 
developed. One particular motivation for doing this is the 
tropical nature of the Mei-Yu front. The intense convection, 
which is correlated with the occurrence of strong low (ieee 
vorticity, points to a different possibility. The large 
latent heat release provides an energy source which may, in 
certain situations, be converted to the kinetic enengy sam 
the LLJ. This possibility is worth investigation especially 
in light of the recent symmetric instability studies, which 
show that slantwise convection may occur in a frontal region 
even if the symmetric instability condition is neutral or 
slightly stable. 

Therefore, the following topics are pursued in this 

study: 

1) Study of the structure and dynamics of a steady 
front, which is investigated by varying the ver tigen 
stability, moisture content, Coriolis parameter and 
boundary-layer fluxes in the model; and 

2) Study of the relationship between the occurrence of 
the LLJ and the convective activity in a Mei-Yu 
tront. In  partveular. the possible role of a 
secondary circulation equatorward of the front caused 


by convection will be explored. 
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In Chapter II, the basic formulation of the model is 
presented. In Chapter III, the results of the numerical 
experiments are described, but the main interpretation of 
the results with respect to the dynamics of the LLJ is 
discussed in Chapter IV. Chapter V gives the summary and 


conclusions. 


Ba 


Fig. 1.1 


Mei-Yu season 





Month 1 3 § ‘3 i: 30 C 31 


Mean daily rainfall (y-axis, mm) for the period of 1951-1970 
at Taichung, Taiwan, which is a typical station having a 
Mei-Yu rainfall maximum (Chen, 1977). 
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Fig. 1.2 Mean onset date of the summer monsoon and the onset date 


of Mei~Yu season in East Asia (Tao and Chen, 1986). 
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Fig. 1.3 Frequency of frontogenesis at 1° longitude by 1° latitude 
grid intervals during the Taiwan Mei-Yu period of May 15 to 
June 15, 1968-1977 (Chen and Chi, 1980). 
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Vertical cross section of (a) east-west wind (u; solid) and 
north-south wind (v; dash) components in m s~}; (b) total 
wind speed (solid) in m s! and mixing ratio (dash) in 
g kg-l along 120°E at 1200 GMT 11 June 1975. J indices 
from left to right increase from south to north with a grid 
distance of 240 km. The Mei-Yu front at 850mb is near J=7 
(Chen, 1978). 
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GMS-1 satellite IR picture at 0000 GMT 28 May 1981. 
and c indicate organized mesoscale convective system. 
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II. NUMERICAL MODEL 


This study uses a modified version of the model 
developed by Williams (1974) and Williams et al. (1981). 
The frontogenesis in the model is forced by a nondivergent 
basic wind field that contains stretching deformation and 
includes an Ekman layer’ structure. The model employs 


vertical and horizontal diffusions of potential temperature, 


moisture and velocity. The diffusion coefficients vary only 
with elevation. Pressure coordinates are used in the 
vertical and a tropopause is simulated. This should lead to 


a more accurate simulation of the upper regions where latent 
heating will be important. The time-dependent quantities in 
the model are functions of y (cross-front) and p and the 
model is integrated until a quasi-steady state is achieved. 
A planetary boundary layer is added with aerodynamic bulk 
formulae for the surface fluxes. A moist convective 
adjustment process ins _ included to represent cumulus 
convection. A staggered grid is used in the horizontal 
because it gives a better treatment of the solutions that 
involve heating. 

In Section A, the Diasane equations in pressure 


coordinates are developed including the appropriate moisture 


terms. The convective adjustment and tatent heating are 
discussed in Section B. In Section C, the boundary layer 


treatment is described. The initial conditions and 


eG 


tropopause simulation are in Section D. The details of the 
specification of model @b lc ceame mere arrangements and finite 


difference scheme are in Appendix A, B and C respectively. 


A. BASIC EQUATIONS 


The hydrostatic primitive equations in pressure 


coordinates may be written 


——<— 


QV 6 . Q OV (2.1.1) 
xy t v-(VV) + 3p (oV) + Vigan ave hegev + ) (K(p)ao J 3 

08 3 a 3 86 (2. ieee 
ae + ¥ (GN) + 3p Gw) = Agv26 + 3p EK(p) ao) Ok eae > 

Q 2 _ 3 38 

se + (QV) + ae (wq) = Aqy*q + 3p [K(p) “rg +Myn - Mc , (2. ae 
VV +dwfAp = 9, ¢ le 
o> /op a -r(p)e, (2 1 5) 


where: 


V = horizontal component of the velocity, 
w = dp/dt, 

p = pressure, 

d = gZ = geopotential, 


6 = T(pp/p)k = potential temperature, 


— 
I) 


temperature, 
R = gas constant, 
Cy = specific heat at a constant pressure, 


K = R/Cy, 
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22 sind = Coriolis parameter, 

rotation rate of earth, 

latitude, 

specific humidity, 

horizontal diffusion coefficient of momentum, 
horizontal diffusion coefficient of heat, 
horizontal diffusion coefficient of humidity, 
vertical diffusion coefficient, 

heating due to dry convective adjustment, 
heating due to moist convective adjustment, 
large scale latent heating, 

Specific humidity change due to moist convective adjustment, 


Specific humidity change due to large-scale condensation, 


h 0 
em Do" P Op \ 7 
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The upper and lower boundary conditions are given by 


w(X, Ys Pts ic) Sean Ys Pos i) = i. 


V(X, ¥s Pos t) = 9, 

ov _ 

pix Ys Pts t) = 0, (2.10 
06 eo 

spk Ys Dt » t) = 5 pte Ys Pos t) = Q, 


Here pt and po represent the pressure at the top and bottom 
Of the atmosphere, respectively, and these quantities are 
approximated as constants. The conditions (2.1.6) include 
the kinematic boundary condition at the top and bottom of 
the atmosphere, the zero slip condition at the surface, the 
zero stress condition at the top of the atmosphere, and the 
zero heat flux conditions at the top and bottom of the 
atmosphere, respectively. 

The approximate steady-state solutions’ to (2.1 2 


(2.1.5) that satisfy the boundary conditions co ieen are as 


follows: 


DIx(1 - eS cos &) + ye sink] 


— 
il 


v = -NLy(1 - eS cos &) + xe% sine] 

w = 0), 

6 =O = D(x? + y)/2 - FDxy 4+F(p), os 
o =e 

dee pae 
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where &= [f/2K(po)]+’2(po - p). The deformation D, and the 
Coriolis parameter f are constant, and the functions F(p), 
G(p) and H(p) are arbitrary except that F and G must satisfy 
mmemayarostatic equatiom (2.1.5). These solutions are only 
approximate within the Ekman layer because some. small 
advection terms have been neglected. 


The dependent variables are subdivided as follows: 


J =e Y5 p) + u(y, Be tiie 


V = v(x, y, p) + v(y, P, t), Gomi) 
wey. pat. © lacy, pret), 


(eo Oem ecm Dp, t),, G =aq(y, p, t), 


where all of the departures from the steady-state solutions 


are assumed to be independent of Stil If we substitute 
meer rewations (2.1.8) into @2.1.1) - (2.1]l.d) and use 
i OV | i au _ ON _ | 

Petes, ~ O(L - eF Gases), amy) “lee 7 ve sinf&, u, v satisfy 
the steady-state equations. The equations of the departure 


field from the steady-state solution become 


Ou BOnL koe: ou, OU 4 ities du 
pt By (uv) + ap (mw) + May Yai y Vay W > 
eee + (om 
May op sein (22. le) 
3S / 


a ae OW 
pet ay (vv) + a (wv) + sy bw) tus os 
(2.1. Ge 
2 0 Ov 
= -fu - "+ An Se a op ‘Kap? 
a0. 3 2 ny 4% 08 = 9 020 5 cy o& 
St oye) ay Oo ae ca ema ap ‘K dp? 
(2.17 
+ yn + Om + He 
2 
OG ipgussORs ace ty oo Aiea” 4-8 (ees 
+ Ma - Me , 
OV OW _ 
sy * dp (2. lean 
On 
— = -f )e 
Op (P (2.1. 


Using the relationship of (2.1.7), the (2.1.9) becomes 


oO 


oe epee noe re 
<F + By (uv) + SD (wu) + D(1 é cos © )u 


- ALy(1 -e cos &) - xe% sin 1x4 Die sin &)v (2.1.15) 


+ od 3 Ex] ~ eS cos ee) + yes Sin & | 
= fy + Ag Oy OK OY) 
MOY Op aici 5 
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foe ismeeto be hoticed that the assumption timat u is 
independent of x for all time is violated in (2.1.15) in the 
two terms where x appears. This indicates that the other 
quantities which were assumed independent of x are actually 
x-dependent, which would Drine other terms Te © the 
equations. In this study, these added terms are neglected 
m@eece!. lo) and the other equations are applied at x-0. If 
the initial u, v, w, 98 and q fields are independent of x, 
the error from this approximation will grow slowly and it 
will be confined to the boundary layer since the x-dependent 
terms in (2.1.15) are zero outside the boundary layer. 
mmumee trontogenesis occurs very rapidly, it is expected that 
the development ue X-variations in the dependent variables 
would only have a small effect on the resulting quasi-steady 
Peont . In any case, these effects could not be observed in 
the atmosphere, because atmospheric fronts always have some 
variation in the basic fields along the front which would be 


much more important. 


When (2.1.15) is evaluated at x = 0, it becomes 
oo on ~ et Ou 
et By (uv) + ap wu) + N(1 - eS cos E)u +P 4 
2 
= sj A ee si : re fy OU 
+ N(esS sinE)v -wVp ap fe sin &) fv + An 5yz * 5p (K 55) 
eon) 


So 


where Vpn = -ly “and aa Vyn(l - eS cose). To remove the 


vertical mean in (2.1.10), the vertical average is defined 
as 
Do 
<( )> = ; — I (. igor 
0 i Dt 


By vertically averaging (2.1.13) and using the boundanys 


conditions (2.1.6), we obtain 


0 
cy > ae 
Oy Y 


This equation states that the total disturbance y mass flux 
ls independent of y. It can be expected that’ at “a Tales 


distance from the frontal zone the disturbance mass” feu 


will vanish. Thus, we get 


Vo =e (2. 1 


If we take the vertical average of (2.1.10) and use (2.1.6) 


and “U2... 7s swer-Gibuauin 


K_Oovy . (2.1.18) 


0 <I> 
= «= ees <u> + Cae ae 
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Wien hugo result 1S subtracted from (2.1.10), we have 


OV 0 0 0 = - 
oe (Sees ~ <yv>) + —— wv) + = (vw - <vv> 
: + sy (w vv>) 5p wv) ae ) 

OV V av OV 

ee on aS ee = > 
po ae <u a naa Wisp 
meee = <y>) eos (x - <>) 

Oy 


Po 
Using (2.1.7) and evaluating at x = 0, we obtain 
dv, d 0 oo = 
cea ay (vv - <vv>) rain (wv) + ay (rv - <v>) 
+ Dues sin& - <ue™ sin &> - Vplu 75 (eS cos &) - wes (eS cos &)>] 
= ale -~ <>) - f(u - <v>) +A o“v 
ay. CO ny 
3 OV K Ov 
+ —— (K —) - (—=— ; Zei.19 
ap MK Sp) - pap), ) 


Poe cel. lly evaluated at x = 0 may be written, 





068 d(ve) 0 06 

=. + en —- 

f° Oy ap OO te ay Cio) 
= 0726 0 06 
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Similarly, (2.1.12) beeemes 


oq, 2 3 3 
att oy VV) + oR a) tT By 


Ae ‘ (2.122 
= Aq Soo + Hp (K(p) S87 + My = Me 


Integrating the hydrostatic equation (2 los) with 


respect to p and removing the vertical mean, we obtain 


p p 
m-<a>=f r(pledp-<f r(p)edp>. (2. lee 
Pt Pt 


With G2 ease) (eZee Ls) and (2s 19) — C2 Slee form a 
complete set of equations which can be solved by a marching 


process. 


B. CONVECTIVE ADJUSTMENT PROCESSES AND LATENT HEATING 
This section describes the dry and moist Conve@isme 


adjustment processes and the latent heat release that are 


used in the model. The dry convective adjustment process is 
activated in any layer where d68/dp > Q. In such a layer, 
8 is set equal to the average value within the layer. -: Then 


the layers above and below are tested for instability and 


modified if necessary. This process is represented by Qa in 
(2.22200 

When the relative humidity exceeds 100%, that is 
q > qs [qs(T,p) is the saturation mixing ratio], latent: 
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of condensation will be released. To adjust the relative 


humidity to 100%, T and q must be changed to satisfy 
q + 5q = qg(T + 5T, p), ee) 
GpbT = -Léq, (2.2.2) 


where L is the latent heat of condensation, and §q and 6§T 
are the changes in q and T during the process. A good 


approximation to the function qs(T,p) is given by 


1,38 E ] l 
ds = —p— exol=a™ apg >); (2.2.3) 


where Rv is the gas constant for water vapor. 

Mite wequations (2.2.1) and (2.2.2) require that the final 
state of the air be exactly saturated and that excess 
moisture be condensed isobarically, which releases the 
latent heat to the air. 

Bauations, (2:2.1) and (2.2.2) can be expressed in 


1terative form as 


6q = de (1 - r) i SI), 67. q = de (2625.4. ) 
and 
CpdT = -Léq . (2.2.5) 


Mmemninatame 90! from (2.2.3) and (2.2.4) gives 
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Ge ei) 
Sq = ” (2 SG 
L tae 
p Pp 


At constant pressure, the Clausius-Clapeyron equation can be 


expressed as 


(FS) ae aS (2 


Substituting (Z. 276) 8tntor(2. Zeer es 


q<e(1 - r) 
5q = {2 : 


et “Eph IE CH 


6q may now be caleulated Metrome 4 2.253) as can the 


corresponding value of 6T from (2.2.4). The first estimates 


—- aio, 


of the adjusted temperature and mixing ratio are, 


T' =T+6T and q' =q+6q. (2). 2p 
These values are in turn used to obtain improved 
approximations. Iterations are continued until the values 


of’ T’ andiaq att each level match with sufficient accuracy a 
value of 100% relative humidity. The changes are then 
applied to the potential temperature and mixing ratio fields 
of the model, which accounts for the term Hc of (2.1.20) 


anGgwMemof (2510 Zit 
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The convective latent heating is parameterized by a 
convective adjustment scheme. The scheme yields similar 
results in describing the large-scale effect of latent 
heating as an explicit scheme which included prediction 
aaueens fOr water vapor, cloud water and rain water (Hsie 
and Anthes, 1984). 

MmciGedom wietemcOnveection 1S OCcUurring, it 1s necessary 
to determine whether the atmosphere is convectively unstable 
to moist processes. A layer will be conditionally unstable 
if 08e/dp > 0, where Be is the equivalent potential 


temperature, which is given by 


om = Cees) (2.2.10) 
Cyl 


As with the dry convective adjustment, @e is set equal to 
the average value in the unstable layer. This gives a value 
for Be at each level and 6T and 6qG can then be 
determined from the following calculations. 

The heat release, Qm, as a result of moist convective 
adjustment is included through knowledge of the change of 
equivalent potential temperature during the adjustment. 


Taking the natural logarithm of (2.2.10) gives 


£nO, =2n6 + os . aa 
p 





A constant pressure surface may be assumed since the model 
uses pressure coordinates. Differentiating both sides of 
(2en2. lela) and recalling the definaticn of potential 


temperature gives 


se 5 del. eo. 
Tf oe -—S— . (2 2ailiae 


Solving (2.2.65) for dqs, substituting into equation (2923 


and solving for dT leads to 


——= 
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1 L L 


where 66e is tne change in @6e needed to adjust the 
gradient d8e/dp to be equal to zero. 

ad atmosphere remains saturated throughout the 
process of moist convective adjustment, the Clausius— 
Clapeyron equation is used again to obtain the corresponding 
6q or actual éqs. Changes are again applied to the 
potential temperature and mixing ratio fields of the model 
which accounts for the terms Qm of equation (2.1.20) and Mc 
of equation (2Z23I2Z2h 2 “The convective adjustment with a 
saturated atmosphere is a rather crude parameterization of 
the cumulus effects. However, the relative humidity can 
indeed become very high in intense Mei-Yu cases (Chen and 


Tsay. (97 Sie Furthermore, the strong control by large-scale 
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frontogenesis forcing makes the results less sensitive to 


the details of the cumulus parameterization. 


C. BOUNDARY LAYER PARAMETERIZATION 

The surface fluxes are parameterized by the  bulk- 
aerodynamic formula and the surface fluxes of heat and 
momentum are assumed to influence directly only the lowest 
model layer. The magnitude of the ice. stress is 
proportional to the square of the wind speed and the surface 
sensible heat flux is proportional to the product of the 
wind speed and the vertical gradient of temperature between 
the ground and the middle of the lowest model layer. ais 
formulation is best suited to models with coarse vertical 
resolution, since the depth of the lowest layer should 
correspond to a typical PBL depth. 

In this study, the vertical sub-grid scale mixing terms 


in Eqs.(2.2.16) are applied to the lowest layer of the 


model. The term may be written as 
Oo Ou fF Ss 0% Ose) 


At the top of the boundary layer, we assume 


Fy (p) = Fy(z) = K(z) 4 ‘ a ee 


K(z) is defined as 


K(z) = Cplul 
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where Co is the drag coefficient and |u| = +¥ ue + v2 ° 


Applying (2.3.2) to the lowest two lbayersy po and prj ot jane 


model, we get 


Fy = -pg Cylul(uz - ug)/(pz - po) 
= - og Colulu,/Ap 
cslmee Uo = er rom the lowes boundary cama ten and Ap] =iiaa - po . 


Therefore, the u momentum flux at the surface becomes 


_ - Ch {ul 
an - aaa hee ) uy, ( Qualls 39 


where 10-3 is a unit conversion factor which converts the 


unit of pressure to cb. Similarly, the v momentum flux is 


Foo 1073 29 Cp jul y,, : 


as Yo" Ap (2.354) 
In this study, the PBL depth, h, is estimated as 
Ay 2 100 - 91.9 
: if . 3 = 
h = ——s4 [oe Oe ee 
The sensible heat and moisture fluxes are 
Fp = - 103 (29MM) @, -0,) 
g§ 7-7 Ap 1 0) (2.35a—) 


= = in (2 fo ) (q, a qy)> (2.3.6) 


4% 


ean — > FF 


where 680 and qo are the specified surface values. 


The latent heat flux is defined as 
; C2 he) 


mmere LF = 2.5 x 105 Joule/kg, cp = 1003 Joule/(°K/kg). 


D. INITIAL CONDITIONS 


The initial potential temperature is given by 
2 ; 
@(y, p, 0) =67(p) - a arctan (sinhay) , ea) 


where 9;(p) is a prescribed vertical potential temperature 


profile and a is one-half the total horizontal potential 


temperature variation. ® is defined as and is evaluated 
, 6 
fx /H —— (—_} (eg) at some middle level p = pn. 
Pm ° Pm Op Pm 


The thermal wind equation, which is derived using the 
hydrostatic equation and the definition of geostrophic wind, 


may be written as 


eae ea. 


meeerprating (2.4.2) from p = pm to p = p 


p p . 
Ou R 06 1; p~p 
J s-dpeese f -(L¥ a 
0 f 0. p 
and assuming u = 0 at p = pm, we obtain 
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o 
Uys p, O) = = SS CB - (Boy (2.4.3) 
0 


USine 802.4) ee 


farctan (sinh ay) |] 


CY 
cD 
I~ 
ala 
> 


Oy 

=. 220 sech (ay) -. 
T 

Substituting (2.4.4) into (2.4.3), the initial wu componmemm 


of the velocity may be written, 


2axC 
u(y, p, 0) = - “Sep at “ae sech (ay) . (2.4.5) 
Q 
The v component of the velocity is set initially to zero, 


1 e e 5 


V-CY. 5 iD uma 


SO jthaiee BiCeg o, O) = 

The initial field of mixing ratio is calculated by 
determining the saturation vapor pressure es at each grid 
point from an integrated form of the Clausius-Clapeyron 


equation. According to Hess (1959), 


L 1 Ik 
eg = 0.611 exp (—- (4=—-- =)) , (2.4.6) 

S oie 
where To 1s a constant of integration. Knowing the 
saturation vapor pressure, the saturation mixing ratio is 
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(2 2a 


a hf lee 


ae — ae — 


as 


obtained from the approximate relation (Haltiner and Martin, 


god). 


Us = es : (e247) 
Knowing the potential temperature at t = 0 from (2.4.1), it 
is possible to calculate qs at each point uSing equations 
m=. 4.6) and (2.4.7). Defining q to be a given fraction of 
qs allows the testing of atmospheres with different initial 
relative humidities. 

The initial conditions described in previous’ sections 
will be constrained by a simulated tropopause that is 
centered at about 300 mb level. The simulation is 
acconplished by using a reference temperature, Trer, at 
about the midpoint of the horizontal scale and calculating 
temperatures such that a tropopause with the expected 
variation in height is determined and a zero vertical 
temperature gradient is achieved in the lower stratosphere. 


A reference temperature, Trer, iS computed using (2.4.1) 


Tref = 215 +a (=) arctan (sinh ay) emeece’)" : COD Tey 
0 | 


When the temperature in each column reaches this reference 
temperature, it is held constant at that value, as expected 
in the lower stratosphere. The temperature is then 


converted to potential temperature. 
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TIi. NUMERICAL EXPERIMENTS AND RESULTS 


In this chapter two groups of numerical experiments are 
carried out to simulate the midlatitude and subtropical 
Situations, respectively. For each’ group, the Coriolis 
parameter and initial temperature sounding are fixed, but 
the amount of the initial relative humidity and surface 
fluxes of momentum, mixing Be and genau heat are 
varied. The purpose here is to describe the structure 
differences due to the difference in latitude, relative 
humidity and surface fluxes. The discussion of the results, 
which focuses on the development of the low level jet, is 
given in Chapter IV. 

We will first give the parameter ranges used in the 
experiments (in Section A), and then describe the numerical 
results of the midlatitude (Section B) and subtropical 
(Section C) cases, respectively. As mentioned in the 
Introduction, the main purpose for all numerical experiments 
is to obtain a quasi-steady solution. However, it is also of 
interest to follow the time evolution of the frontogenesis 
process to see how the quasi-steady states are developed. 
Therefore, in each group the evolution for a "typical case" 
1s shown. It turns out that these oSoimthrone provide useful 
insights about the different nature of the effect of cumulus 


convection in the midlatitude and subtropical cases. 
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A. EXPERIMENT CASES AND PARAMETERS 


The following values are used for the constants in this 


study: 


9.81 ms72. 

1903 joule/kg/°K 
287 joules/kg/°K 
R/Cp 

461 joules/kg/°K 


2.5 x 10° joules/kg 


5°K/km 


- 0.53448°K/cb 


100 ch 


55° Cb 


Ag = My = 
Kq = Ko 


74) Sec 


3 x 10" WW /s 


0.0007 ch2 /s 


60 km 
2760 km 


Each 


experiment is designated by a three-letter case 


name as follows: 


LRF 


where L defines the latitude environment. 


ees 


M refers to a midlatitude environment, with the 


following parameter settings: 


a. 


b. 


Zeal 


f = 1.06 x 10-4 s-! (corresponds to 46s 
initial ween cal potential temperatiiies 
distribution given by a typical midlatiGmae 
sounding (Fig. 3.1); 

initial meridional temperature variation 1s 
given by a = 8°K; and 

large-scale deformation fields given by D = 
lx TOs oss * . 


S refers to a subtropical environment, with the 


following parameter settings: 


a. 


De 


f = 0.558 x 107-4 s~! (corresponds to 22.,5°%e 
initial vertical potential temperature 
distribution given by a typical subtropical 
SOUNndI WoC! 1 ooo le 

initial meridional temperature variation is 


given by a = 4°K; and 


large-scale deformation fields given by D 


0.5 x GOR woes 
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3. R defines the relative humidity: 


a. R= 6 means R.H.=60%; and 


i 
tl 


Bietueeeealt 9 means R.H.=90%. 

The observed low-level relative humidity in the vicinity 
of a Mei-Yu front varies between 60% - 90% (Chen and Tsay, 
1978). When the Mei-Yu front is north of the Yantze River (a 
midlatitude situation) the minimum R.H. is about 60%. When 
it is located over Taiwan — Sesneasweun coast of the 
Chinese mainland the maximum R.H. is near 90%. 

4. F defines the surface fluxes: 

a. F = 0 means there is no surface fluxes of heat, 
moisture, and momentum; and 

b. F = 1 means surface fluxes of heat, moisture and 
momentum are included. 

Typically, the fluxes are large when a subtropical 
portion of the front is over a warm sea surface. When the 
front is moved over land in the midlatitude cases, the 
fluxes are normally weaker with the exception of momentum 
fluxes. For self-consistency, the fluxes of all quantities 
will either be included or excluded in the numerical 
experiments. For all cases, the initial zonal wind is 
calculated from the potential temperature distribution using 


the thermal wind, and the initial v component of the 


velocity is set to Zero. 
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B. MIDLATITUDE SIMULATIONS 


1. Typical Midlatitude Case (M60) 





The midlatitude case of M60 with an initial R.H.=60% 
and no surface fluxes describes the development of the front 
in detail. The potential temperature and mixing ratio at 
t=0 are shown in Fig. 3.2a, and the zonal velocity at t=0 is 
shown in Fig. 3.2b. 

Fig. 3.3 depicts the deve lomment of the potenti 
temperature and mixing ratio fields. The horizontal 
gradient of both fields increases steadily through the first 
two to three days asa frontal structure develops near 42°- 
43°. The maximum surface temperature gradient nearly 
triples (from the 2°K/degree latitude to about 6°K/degree 
latitude), while the maximum surface mixing ratio gradient 
increases nearly ten-fold (from about 0.5 g/k¢/ deenaee 
latitude to about 5 g/kg/degree latitude). This midlatitude 
mixing ratio gradient is about twice as strong as that 
observed in a subtropical Mei-Yu front by Chen (1978) as 
shown in Fig. 1]1.4b, and the average mixing ratio is lower 
(8 g/kg versus ll g/kg). The axis of maximum mixing ratio 
exhibits a conspicuous poleward vertical tilt from a at 
the surface to 48° at 8 km. This poleward vertical tilt 
Reet a slope of approximately 1.1:100 agrees well with the 
typical midlatitude fronts (Byers, 1974). The tilt is 
slightly stronger than that observed in the subtropical Mei- 
Yu front (slope 1.3:100) aS shown an” Frese tse The 
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development of the frontal thermal structure, including the 
concentration of the isentropes in the lowest 2 km and the 


poleward tilt of the frontal zone from the surface to the 


upper troposphere, is clearly indicated in Fig. 3.4. In 
addition, a secondary maximum occurs in the upper 
troposphere (6 - 8 km or 600 - 400 mb) near 50°, and a 


reversed tilt and reversed gradient is present in the lower 
stratosphere. This — structure is a result ae 
the imposed condition of a peop ause and is consistent with 
previous results obtained by Williams et al. (1981). 

The vertical velocity field (Fig. 3.5) quickly 
develops and becomes organized, with a concentrated upward 
Pmo@tion of > -4.0 x 1074 cb/sec occurring at day 2 in the 
lower troposphere near 43°. A poleward vertical tilt of the 
area of upward motion develops and becomes more obvious with 
time. After day 4, a quasi-steady state clearly delineates 
both the upright maximum rising motion centered at 43° to 
44° in the lower-middle troposphere and the vertical tilt, 
which is well represented by the w=0 line poleward of the 
maximum rising motion. 

A very small area with a secondary maximum gia 
vertical velocity appears at day 2 near the surface ahead of 
ChemPrent (41°). This represents the effect of the forcing 
by the Ekman layer structure of the mean flow, as was 
discussed by Williams et al. (1981). Because this area is 


very close to the surface, it does not cause condensation 
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and latent heat release even in the highly moist cases 
discussed below. 

The initial zonal wind (Fig. 3.2b) is in geostrommae 
balance with the temperature field, with an upper 
tropospheric westerly ‘maximum of ~55 m/s and a surface 
easterly maximum of -25 m/s. The pattern is almost 
symmetric with respect to 46° latitude. As the integration 
proceeds (Fig. 3.6), Aswan develop and a strong oNGeie 
begins to appear at day 1.0 below 6 km. The shear quickly 
reaches a steady state with a frontal structure which tilts 
from 42° at the surface towards 45° at 4 km. This is in the 
same vicinity as the Ww =O line and the frontal zone that is 
defined by the potential temperature and. mixing ratio 
gradients. At the latitude of the surface front, a. local 
westerly maximum develops near z=4 kn. By day 2 “Give: 
3.6b), this narrow westerly regime penetrates downward and 
appears to be an intrusion into the lower tropospheric 
easterlies. It remains in the 3-4 km layer — day 4. 
Although the magnitude of westerly wind in this narrow zone 
is small, the poleward meridional shear is large, which 
results in a vorticity in this region that is comparable oO 
the values near the upper tropospheric westerly jet. The 
development of this narrow westerly regime is a result of 
the intensifying meridional temperature gradient which is 
forced by the frontogenesis process in the lower 


troposphere. Development of this zonal velocity field is 
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also dramatically illustrated by the evolution of the 
relative vorticity (Fig. 3.7). The initial flow pattern 
develops rapidly into a quasi-steady state frontal structure 
after about two days. Maximum cyclonic vorticity in the 


lower troposphere is concentrated on the warm side of the 


Prone. To illustrate the vertical structure of the 
horizontal wind shear with that of the horizontal 
temperature gradient, 020/ d * =m /1OGmKn “for day 4 ne 
merotted im the vorticity figure. A coherent vertical tilt 


of the frontal structure is obvious. 
2. Other Midlatitude Cases 

The effects of varying the initial relative humidity 
(M90), and the inclusion of surface momentum, heat and 
moisture fluxes (M61 and M91) are illustrated by the 
meemiting strueture at day 4 (Figs. 3-°8-3.12). In addition, 
the precipitation and heating rates are also shown in Figs. 
owes and@3.h, respectively. 

a. Effect of Moisture 

(1) No Surface Flux Cases (M60 vs. M9O) 
It is clear that the initial R.H.=90% case 

(Fig. 3.8b) has higher mixing ratios in general and a warmer 
upper troposphere between 40°-50° (best represented by 
comparing the 98=350°K line). The higher temperature is a 
result of a larger latent heat release in the frontal region 
due to the increased moisture content, as can be seen by 


comparing the precipitation rates oP ioe Sie lesa The 
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meridional potential temperature gradient in the surface 
front for the R.H.=90% case is also stronger (Figs. 3.94 
versus 3.9b). The frontal zones as defined by the region 
9@/dy < -0.02°K km-! are quite similiar in width and 
vertical tilt. Above the surface front position (~42°), the 
mixing ratio in the R.H.=90% case (Fig. 3.8b) also shows a 
significant small scale structure with a wavelength near 2 y 
in the 2-7 km layer. This structs is —" less evident ie 
the R.H.=60% case (Fig. 3.8a). This small scale structure 
is near the region of maximum rising motion above the 
surface front position. The stronger vertical velocity for 
the R.H.=90% case (maximum rising motion > 6x10-4 cb/s 
versus 4xl0-4 cb/s for the R.H.=60% case) is consistent with 
the stronger cumulus convection. This convection arises 
because of the existence of conditional instability which is 
due to the prescribed 90% relative humidity. The cumulus 
convection activity appears EG cause the small-scale 
variations in the mixing ratio field. 

In addition to a larger maximun, the 
vertical velocity of the R.H.=90% (Fig. 3.10b) case also 
displays a more complex structure, with two concentrated 
regions of rising motion. The primary region near 42°-43° 
extends throughout most of the depth of the troposphere from 
about 1] km to 11 km. The secondary region occurs in the 
middle-upper troposphere (5-11 km) just poleward of the 


primary region at about 44°, In addition, a third, 
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relatively weak, maximum upward motion occurs’ farther 
poleward near 47°. This structure is referred to as a 
multi-band, slanted rising motion in later discussions. This 
multi-~band structure in the vertical velocity may be 
compared to some of the observed multi-band structure of 
cloud and precipitation associated with cold fronts during 
mre period of CYCLES (CYCLonic Extratropical Storms) Project 
(Parsons and Hobbs, 1983). | | 

As may be expected, the stronger warm core 
and stronger convection of the R.H.=90% case is related to a 
larger latent heat release. In Fig. 3.14, the heating 
distribution generally, resembles the rising motion field 
(Fig. 3.10). This is consistent with a near balance between 
adiabatic cooling and diabatic heating. The R.H.=90% case 
(Fig. 3.14b) has a maximum heating rate of >6 °K/day, which 
is much higher than the 1] °K/day rate of the R.H.=60% case 
( Paro eee This six-fold increase in heating ‘ies 
consistent with the increase in precipitation rate shown in 
Fig. 3.13. The greater production of kinetic energy implied 
by the stronger rising motion of warm air gives rise to a 
more intense horizontal circulation (Figs. 3.1llb versus 
Seeeba). The upper tropospheric westerly jet (Fig. 3.11b) 
reaches a maximum of > 50 m/s (versus the > 35 m/s for the 
R.H.=60% case) near 52°, which is two degrees poleward of 
the jet core of the R.H.=60% case. A small region of weak 


easterlies that develops equatorward of the westerly jet 1s 
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not observed in the R.H.=60% case. It is to be reqcaliied 
here that these are perturbation wind fields; the basic flow 
is westerly throughout the upper troposphere. The 
distribution of the perturbation, zonal wind iS,.conS Tiga 
with an anticyclonic circulation above the warm core that is 
due to heating. The existence of this easterly regime 
limits the middle-upper tropospheric westerlies to poleward 


of the surface front position, and makes the middle 


tropospheric westerly regime appear as a narrow jet-like 
structure that penetrates downward to below z=3 km. The 


resultant low level westerly regime is sharper and more 
extensive than that of the R.H.=60% case. 

On the other hand, the strength of the 
surface easterlies of the M60 and M90 cases is Very Sigg 
This is probably because most of the impact of a higher 
relative humidity is in deep cumulus convection, with the 
resultant increased latent heat release occurring mainly in 
the upper and middle troposphere, as shown in Fig. 3.14. 

The maximum Cyclon Cc vort Loa 1s 
concentrated within the frontal zone in both the upper and 
the lower troposphere for both cases (Fig. 3.1l2a versus Fig. 
3.12b). The stronger horizontal circulation in the upper 
troposphere for the R.H.=90% case is also evident in the 


Wore leony f avec 
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(2) Surface Flux Cases (M61 versus M91) 

When surface fluxes of momentum, moisture 
and heat are present, the differences between the _— moist 
(R.H.=90%) and the drier (R.H.=60%) cases are basically 
Similar to those between the no-flux cases. The R.H.=90% 
case has a warmer upper tropospheric core, as represented by 
the decrease of height of the 69=350°K line, and again has 
the small-scale variation in eaethc eee in the piste 
motion region over the surface front location (Fig. 3.8d 
versus Fig. 3.8c). The meridional potential temperature 
gradient is stronger for case M91 than for case M6] (Fig. 
3.9d versus Fig. 30'c )* The maximum upward motion, 
precipitation and heating rate are much larger for M91 than 
tom «M651 (Fi goes . Ow 3.13 and 3.14). Stronger’ wind 
circulations (Fig. 3.lld versus Fig. 3.llc) are shown in the 
R.H.=90% case. The circulation includes a more intense 
westerly regime of > 50 m s7! in the upper troposphere, the 
development of a narrow easterly regime equatorward of it, 
and a deeper downward extension of the upper westerlies into 
the low-level easterly jet. The differences in the 
magnitudes of vorticity (Fig. 3.12d versus Fig. 3.12c) are 
similar to other fields with much stronger evelonic 
vorticity within frontal zone for the R.H.=90% case. 

It is to be noted that the small-scale 
variations in all fields for the case with high mixing ratio 


and surface fluxes (Fig. 3.8d-3.12d) exist near surface to 
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the north of {renin As %@an be seen an@Fig. 3.15, ‘theme is 
the cold air region where the sensible and latent heat 
fluxes are large and apparently cause a See layer of 
conditional instability. 

On the other hand, there is one noticeable 
difference between the comparisons of the with-flux cases 
and those of the no-flux cases. The multi-band vertical 
velocity structure in the Oe (R.H.=90%) case (Fig. 
3.10b) is not as conspicuous in the with-flux (R.nv=Ssiee 
case (Fig. 3.10d). However, a weak secondary maximum region 
does appear in the latter near 7 km at 46°, which indicates 
a double-maximum structure and a poleward tilt in the 
vertical structure. The physical implication of Sms 
slantwise structure in the upward motion field will be 
discussed later in Chapter 4. | 

b. Effect of Surface Fluxes 
(1) R.H.=60% Cases (M60 vs. M61) 

The inclusion of surface fluxes (Figs. 3.8a 
and 3.8c) does not lead to any discernible differences in 
potential temperature and mixing ratio distributions, except 
the mixing ratio isolines at the surface have been spread 
towards cold air when’ surface fluxes are included. 
Comparison between Figs. 3.9a and 3.9c shows that the 
surface meridional potential temperature gradient for the 
case with fluxes substantially weakens, although the general 


pattern delineating the frontal zone remains the same. 
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These reductions in the potential temperature and mixing 
ratio gradients at the surface are due to the distribution 
of the surface fluxes given by the bulk formula, which 
provides larger heat and moisture fluxes poleward of the 
surface front in the cooler and drier region (Fig. 3.15). 
The inclusion of surface fluxes causes a 
slight increase in the vertical velocity (Fig. 3.10a and 
a LOc:,) . This is expected sine more latent heat release 
occurs with the increased supply of heat and moisture from 
below, as shown in the comparison between Fig. 3.14c and 
Hag. 3.l4a. Since the resultant upper tropospheric 
temperature is not higher (Figs. 3.8a and 383.8c), the 
additional latent heat release in the area of rising motion 
iemecompengated almost entirely by the increased adiabatic 
cooling associated with the stronger upward motion. 
Increased conversion of available potential 
energy, as implied by the in-phase relationship between the 
vertical velocity and potential temperature fields, leads to 
a stronger upper tropospheric jet (40 m/s versus 35 m/s in 
peemeno-ftlux case) (Figs. 3.lla and 3.llc). On the other 
hand, a decrease of the surface easterlies from 40 m/s to 
30 m/s occurs due primarily to the dissipation by the 
momentum flux transfer. The maximum easterly jet core is 
raised to a level slightly above the surface, and the 


strength and extent of the narrow lower’ tropospheric 





westerly region appear to be slightly reduced for no-flux 
case. 

Near the surface, .the cyclonic” vort emia 
decreases in the with-flux cases (Figs. 3.1l2a and 3. T2age 
This indicates that the weakening of the easterly maximum 
and thermal gradient of the front at the surface does lead 
to a corresponding reduction in the horizontal wind shear in 
the lower troposphere. 

(2) R.H.=90% Cases (M90 vs. M91) 

Similar results due to inclusion of the 
surface fluxes are obtained in the comparison between the 
with-flux and no-flux cases for R.H.=90% as in the R.H.=60% 
cases: 

1) The difference in the distribution of potential 
temperature and mixing ratio is not discernible, 
except near the surface to the north of front (Figs. 
3. Sbaianiclies esi 

2) At the surface, the inclusion of moisture and heat 
fluxes weakens the meridional gradients of mixing 
ratio (Figs. 3.8b and 3.8d) and potential 
temperature in the frontal zone (Figs. 3.9b and 
3.9d) because of the stronger fluxes in the cooler 
and drier regions; 


3) The surface momentum transfer leads to a weakened 
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SunbPacersecoctiermlimpomaxamumes(Pags.  3«61llb- -and »8.l1l1d), 
with the easterly jet core occurring slightly above 
the Sa aoe 

4) The inclusion of surface fluxes causes a stronger 
and more concentrated vertical velocity (Fig. 3.10) 
and heating rate (Fig.3.14) in the frontal zone, so 
that the upper tropospheric westerly wind, the 
westerly regime in the =rawale ante lower @MopospWere 
above the surface front are stronger, and the 
meridional wind shear in the lower troposphere is 
weaker (Figs. 3.llb and 3.1lld and 3.12b and 3.1]2d); 
and 

Oye For RoH2=36%, the multi-band structure in vertical 
Velocity (iia 37 0beand 3.10d), for the with-flux 
case is less aeeeuous than that for the no-flux 


case, as was noted earlier. 


Cc. SUBTROPICAL SIMULATIONS 
1. Typical Subtropical Cases (S90) 

To examine the time evolution of subtropical 
experiments, the case of S90 (initial R.H.=90% and no 
surface fluxes) is chosen as the "typical subtropical case”. 
The initial potential temperature and mixing ratio are 
depicted in Fig. 3.16a and the zonal velocity is shown in 
Poe 3. Ob. The choice of the typical subtropical case is 


motivated by an intention to represent the major differences 
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between the mean midlatitude atmosphere and the  wmean 
subtropical atmosphere during the Mei-Yu season. The mean 
ener Opnierl Mei-Yu environment is conditionally unstable 
with a R.H.  ~aroumd 65%. Furthermore, during intense 
convection periods when a LLJ is usually observed, the R.H. 
is around 90% or higher (Chen and Tsay, 1978). While these 


intense episodes are local and transient, and therefore 


cannot be properly studied by the present two-dimensional 


model, they are the only periods when the LLJ is clearly 
defined. Therefore, the somewhat higher value of R.H.=90% 
is chosen so that its results may also shed light to the 
situations when an intense LLJ is developed. In this view 
the "typical subtropical case" may also be considered to 
represent the "typical intense LLJ case." 

As in the eel midlatitude case, the large-scale 
deformation forcing causes an accumulation of the potential 
temperature and moisture isolines to result in a frontal 
structure. However, the approach towards a quasi-steady 
state is slower. Furthermore, all fields develop 


significant small-scale variations on the equator side of 


the front. (For example, see the potential temperature and 
mixing ratio fields shown in Figs. 3.lfasabor These 
variations, which are similar to those in the _ strong 


convection cases of midlatitude experiments, are apparently 
a consequence of intense cumulus convection in the model 


resulting from the very moist and unstable atmosphere. 
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It should be noted that once the initial vertical 
Ge profile is set from the temperature sounding and 
relative humidity, it will remain fixed in time at Laeere1 
boundaries. Since a conditionally unstable profile occurs 
in the lower latitude regions (because of the initial 
potential temperature and mixing ratio distribution), the 
instability will persist on the equatorward boundary 
regardless of convective Brae Ulett 1 Sant in ane interior. This 
slower development into a quasi-steady state is probably due 
to two factors: the weaker large-scale deformation forcing 
(D=0.5x10-5s-1) prescribed for the subtropical cases and the 
presence of strong conditional instability which implies 
more prominent cumulus convective activity. | 
The small-scale features also dominate the vertical 
velocity field such that the frontal structure is obscured 
by high frequency oscillations. To remove these small-scale 
oscillations, all the subtropical experiments fields are 
subjected to a daily average over 32 outputs at 45 minute 
intervals. Two other alternatives were considered: to 
increase the control of large-scale forcing, or to decrease 
the initial relative humidity. However, both approaches 
will reduce the contrasts between the midlatitude and the 
subtropical cases especially those cases with clearly 
defined LLJ, which is one of the main objectives of this 
study. Therefore, it was decided that the differences 


between the midlatitude and subtropical conditions should be 
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kept as large as possible to identify the underlying 
dynamics of the development of the LLJ. Intermediate cases 
with R.H.=70% and 80% will also be run to ee nieie a 
physical significance of the intensively convective R.H.=90% 
cases. 

The daily averages of the potential temperature and 
Mixing ratio, and the meridional potential temperature 
gradient for days 2-5 are Sen in Figs. 3.18 and 2 ore 
respectively. Here each "day" represents the average of the 
previous 24 hours (i.e., day 5 is the average between 96- 
iZ0Uhoe In these diagrams, the small-scale features have 
been filtered out, and the development of the frontal 
structure is clearly evident. The surface mixing ratio 
gradient and range agree reasonably well with Chen’s (1978) 
observations (Fig. 1.4b). 

The meridional potential temperature gradient also 
evolves into a frontal structure near the surface around 16° 
after 4 days (Fig. 3.19). The tilt of the frontal zone is 
Similar to the midlatitude case (Fig. 3.4), but the well- 
defined frontal zone is confined to the lower troposphere. 
Chen and Chang’s (1980) comparison between different 
segments of an elongated Mei~-Yu front that extends from the 
subtropical southern China to the midlatitude central Japan. 
They observed that the subtropical segment of the front has 
less vertical tilt and is shallower than the midlatitude 


segment. These differences were attributed to the stronger 
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Gemvection in the subtropics. It will be shown later that 
in the subtropical with-flux case (S91) the tilt is indeed 
less while the vertical extent remains shallow. 

By day 2, the concentrated rising motion (Fig. 3.20) 
develops in the equatorial latitudes and then moves poleward 
towards the area of maximum deformation forcing. At day 3, 
the maximum rising motion is at 16° with a magnitude 
~17x10-4 cb/sec. After day 4, it shifts to 18° and the 
maximum decreases to 4-5 x 1074 cb/sec. This is the process 
for the model to develop the Mei-Yu front from an initial 
zero-perturbation- state. In the real atmosphere, the 
process is much slower, where the monsoon trough, which is 
developed over a seasonal time scale, is already in place 
when the front moves into the Mei-Yu position from the 
Mert . The development of a slanted structure in the — 
(z=0-3 km) and upper parts (5-10 km) of the concentrated 
rising region is evident. The slanting is northward with 
height, in the same sense as the slope of the frontal 
surface. The multi-band structure is not clearly seen in 
the 24-h averaged vertical velocity field but is more 
clearly visible when instantaneous (not daily-averaged) 
vertical velocity fields are examined. An example is given 
in Fig. 3.21, which is the instantaneous vertical velocity 
at the time-step of day 4.5. 

The evolution of the zonal wind is illustrated in 


Eee S522. The gross structure of an upper tropospheric 
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westerly jet, poleward of a lower tropospheric easterly jet 


with a maximum at the surface, roughly resembles that of the 


typical midlatitude case (M60, Fig. 3.6). However, there 
are significant differences. The subtropical westerly jet 
core is lower (at 10 km), shallower and wider in the 
meridional direction. These features are different from 


what would be expected for the mean polar front jet and 
subtropical jet. Compared to Chen’ é (1978) observation of P 
subtropical Mei-Yu case (Fig. 1.4a), the jet core is lower. 
The reason for the lower jet core (at 10 km versus 12 km in 
fig. 4a) iS apparently the fixed tropopause which is 
imposed by the mean sounding at 10 km. In an intense Mei-Yu 
period, the deep convection and the tropopause are a couple 
of kilometers higher, where the jet maximum is also located. 
On the ener hand, the meridional extent of the simulated 
westerly regime and that observed are comparable. A 
significant easterly regime occupies the entire tropical 
upper troposphere and caps a westerly regime in the tropical 
lower troposphere. In Fig. 1.4a, the observed zonal wind is 
also easterly in the upper level tropics, although the area 
is somewhat higher and more equatorward than the simulation. 
The imposed tropopause in the model is again the reason for 
this difference. 

Within the tropical low-level westerly regime is a 
remarkably strong and narrow jet structure centered at 2z=3 


km at the latitude of maximum rising motion. This low level 


Ze 


lett leet appears atveedaye 2 withea magnitude ~5 m/s. et 
develops to > 10 m/s at day 5. In terms of the relative 
position with oe to the Bit this low level jet is 
Similar to the local lower-tropospheric westerly regime 
simulated in the midlatitude cases, but the intensity is 
much stronger. This intensity agrees reasonably well with 
the observed LLJ strength in Fig. 1.4b, especially when the 
addition of a mean westerly is BRSERRIOY cae | 

The strong horizontal wind shear that develops 
anoundethe frontal wones is»shown in Fig. 3.23. Compared to 
laeg . Biel. this subtropical case has a deep lower 
tropospheric shear zone that extends into the middle 
fe Sphiener and its magnitude poleward of the low level jet 
Significantly exceeds that of the upper tropospheric 
maximum. The slope of the frontal shear zone is again 
larger (less vertical tilt) in the subtropical case than in 
the midlatitude case (Fig. 3.7), at least in the lower 
troposphere. In Fig. 3.23d, a region of strong negative 


relative vorticity exists above the surface front near 


10 kn, which indicates the possibility of inertial 
inswab teh t ye Des 6} see thesesconsequencesm of the intense 
convection and may not be physically realistic. It is not 


known whether this instability contributes to the small- 
scale variations. 
To demonstrate that the five-day integration gives a 


reasonable quasi-steady state, the averages of each field 


we 


over the first five days of integration was used as the 
initial conditions for another integration of 120 the The 
results of acuySee ven temperature and mixing Yabaor 
meridional potential temperature gradient, vertical velocity 
and zonal wind for the (daily averaged) new day 5 are shown 
in Fig. 3.24a—d. These fields are similar to those of the 
original daily averaged at day 5 (Figs. 3.18, 3.19, 3.20 and 

Seeei2 ).2 This result indicates ana the daily averaged rae 
are stable in time, in spite of the presence of small-scale 
variations and possible inertial instability. 

The subtropical frontal positions as indicated by 
the various fields are farther equatorward from the central 
eae of deformation forcing than that of the midlatitude 
cases. In the midlatitude case, the central axis of 
deformation is at 46° and the frontal structure near the 
surface develops at about 42°. In the subtropical case, the 
axis of deformation is at 22.5°, while the surface frontal 
zone is near 16°. This increased distance ap pale 
reflects the effect of the smaller Coriolis parameter which 
leads to a relatively larger ageostrophic secondary 
‘circulation around the frontal zone. The advection” byaiaa. 
lower branch of this enhanced secondary circulation tends to 
move the surface front farther equatorward. In adden 
the effect of a possible secondary circulation emai of the 


front, that is induced by the strong convection, may also 
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comer bwte*™ tor the more equatorward position of the surface 
ierroult. 
Zeeoruner Subtropical “@ases 
As in the midlatitude experiments, the effects of 
varying the initial relative humidity and surface fluxes are 
studied by running three additional subtropical cases (S60, 
molt and S91). The resultant structure at day 5, averaged 
over the preceding 24-h pari ge for all i ovr subtree cal 
Gases are shown in Figs. 3.25-3.29. The precipitation and 
heating rates are shown in Figs. 3.30-3.31. In these cases, 
no significant small-scale variations are present, so the 
averaging is unnecessary. Nevertheless, it is done to make 
the” analysis procedure of all subtropical cases consistent. 
a. Effects of Moisture 
Panels a and b in Figs. 3.25-3.29 are for S60 
CReH.=60%, no-flux) and S90 (R.H.=90%, no-flux) cases, 
respectively. The R.H.=60% case (Fig. 3.25a) has a less 
intense moisture front throughout the lower troposphere. 
The poleward vertical tilt of the maximum mixing ratio 
gradient is about the same as in the R.H.=90% case (Fig. 
ome25b). However, the curvature near the axis of the 
maximum mixing ratio does not appear as sharp as in the 
R.H.=90% case due to the drier air and weaker gradient. The 
short- wave feature in the middle troposphere over. the 
surface front position is not as intense as in the R.H.=90% 


case, but it is still very notable, which indicates the 


digs. 


presence of Significant cumulus “eonyeeiamen- As may be 
expected, the R.H.=60% case has a lower temperature in the 
upper troposphere (Figs. 3.25a and 3.25b). The small-scale 
variations near the surface in the cold air (Fig. 3.2008 
3.27d), as seen in midlatitude cases, evidence the shallow 
layer of conditional instability caused by surface heat 


fluxes. The conditional instability can be examined further 


by calculating the Oe field which will beseshownedatemmenmn 


Fig. 4.6 (midlatitude cases) and 4.7 (subtropical cases). 
The weaker meridional potential temperature gradient in the 
frontal zone can be seen by comparing Figs. 3.26a and 3.26b, 
both of which also show less small-scale features elsewhere 
in the domain. The vertical tilt of the frontal zone is 
somewhat larger and the frontal zone extends to a slightly 
lower height at the poleward boundary. 

The vertical velocity of the R.H.=60% case (Fig. 
3.27a) has a concentrated zone of rising motion in the lower 
and middle troposphere, but the magnitude is less than in 
the R.H.=90% case. 

The zonal wind of the R.H.=60% case is much 
weaker than in the R.H.=90% case (Fig. 3.28a versus Fig. 
3.28b). The upper tropospheric westerly jet maximum of 
< 25 m/s is only about 70% of the > 35 m/s magnitude of the 
R.H.=90% case. This is consistent with the expected 
difference in the energy conversion between the two cases, 


as the heating rate and the upward motion (which correlates 
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almost completely with the heatin)g of the warm air is 
stronger for the R.H.=90% case. The former is also situated 
more equatorward. The low-level westerly jet structure 
still persists, but the magnitude is reduced by half. The 
upper tropospheric easterly regime almost disappears, and 
the surface easterly maximum is only reduced slightly (Fig. 
ewes and Fig. 3.28b). These features also appear in the 
comparison between Figs. 3.29a — by, which pea 1 y 
indicates a significant reduction of upper tropospheric 
horizontal wind shear. 

When fluxes are included, the comparison between 
the R.H.=60% (S61) and R.H.=90% cases (S91) [panels c and d 
in Figs. 3.25-3.31, respectively] yields similar results. 
In general, the meridional temperature and mixing ratio 
gradients and the zonal wind are stronger in the R.H.=90% 
‘case. —_ 

The subtropical R.H.=90% cases are included to 
study the intensely convective Mei-Yu situations when the 
LLJ is fully developed. However, the high initial moisture 
and the subtropical sounding make the atmosphere highly 
conditionally unstable. The resultant high frequency, small 
scale variations are removed by the 24 h averaging. Even 
though the averages remain stable for extended integrations, 
a question may still BG asked as to whether the results are 
not distorted seriously by the small-scale variations. 


Therefore, two intermediate with-flux cases, R.H.=70% (S71) 
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and R.H.=80% (S81), are run to assess’7 the physical 


significance of the S9l case. Fig. 3.33 compares the 24-h 


averaged zonal wind of all four subtropical, with-flux 
cases. Here it can be seen that the pattern of all four 
cases are similar. The midlatitude upper level westerlies, 
lower level easterlies, tropical upper easterlies are 


present in all cases, and the intensity of these features 
increases systematically as ne initial Re increases Piaf 
60% to 90%. The LLJ is indicated by a very small area of 
5 m/s in the S6l case (Fig. 3.33a), this area expands 
considerably with increasing humidity. When the R.H. 
reaches 90% (Fig. 3.33d), the jet maximum becomes > 10 m/s. 
The similarity in the structure, and the consistent 
intensification of the circulation with increasing moisture, 
indicate that the subtropical R.H.=-90% results can indeed be 
viewed as physically relevant to the strong LLJ situations. 
b. Effects of Surface Fluxes 

When surface momentum, heat and moisture fluxes 
are included (Fig. 3.32), these important contrasts continue 
to exist, as seen by comparing panels b and din Figs. 3.25- 
Bes Ly. Similar to the midlatitude experiments, the addition 
of surface heat, moisture and momentum fluxes leads to a 
higher upper tropospheric temperature. This is apps tae 
due to the ic Hoaeen precipitation and the associated latent 
heat release (Figs. 3.30 and 3.31), which are to be expected 


from the increased sensible and latent heat supply (Fig. 
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3.32). A weaker potential temperature and mixing ratio 
gradient at the surface (Figs. 3.25 and 3.26) is due to the 
stronger surface fluxes of heat and moisture in the cooler 
and drier region (Figs. 3.32). 

The vertical velocity of the with-flux case 
displays the stacked multi-band structure slanted upward and 
poleward in the surface frontal zone (Figs. 3.27b and d). 
Also, the Bini eu of apie is greater than in the wifes 
flux case. This is consistent with the precipitation and 
heating rates (Figs. 3.30 and 3.31) and is due to the 
additional moisture supply from the lower boundary (Fig. 
he 3 Mh. 

Comparisons between Figs. 3.28a and 5. 2OC wand 
between Figs. 3.28b and 3.28d show the effects of surface 
fluxes on the horizontal circulation, which are basically 
Similar to the results of the midlatitude experiments. When 
fluxes are included, the upper tropospheric zonal winds, 
including both the westerly jet and the less intense 
easterly jet, are enhanced. The lower tropospheric westerly 
jet near 3-4 km for the R.H.=90% case (Fig. 3.28d) also 
intensifies. Due to the momentum dissipation at the surface 
(Fig. 3.32), the surface easterlies weaken, and the jet core 
is shifted slightly upward (Figs. 3.29b and d). 

The comparisons of the surface fluxes and no 
fluxes for the R.H.=60% cases (panels a and c of Figs. 3.25- 


3.29) basically illustrate similar results. 
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SUMMARY OF THE RESULTS 


The results of the foregoing numerical experiments may 


be summarized as follows: 


1) 


2) 


3) 


In the midlatitude case, the zonal velocity field is 
dominated by a westerly regime occupying almost the 
entire upper troposphere and an easterly regime 
occupying the lower ropes pee On the other hand, a 
the subtropical cases, the zonal flow equatorward of the 
front tends to reverse, with expansive easterlies in the 
upper troposphere and weak westerlies in the lower 
troposphere. A small regime of easterlies develops in 
the” “madlat i (ude spc ic troposphere only in the high- 
humidity cases. 

At the latitude of the maximum rising motion, a narrow 
westerly wind regime occurs in the lower troposphere. 
In midlatitude cases, it appears as a small downward 
(and equatorward) penetration of the upper tropospheric 
westerlies. In the subtropical cases, it exists 
distinctly as a well-defined, westerly maximum similar 
to the low-level jet. 

The poleward tilts of all fields except vertical 
velocities are more prominent in the midlatitude cases 


are either comparable or stronger than in the 
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5 ) 


6 ) 


subtropical cases. The frontal zone also extends higher 
in the midlatitude cases. These are in agreement with 
observations (e.g., Chen and Chang, 1980). 

The inclusion of higher initial humidity values or of 
surface fluxes tends to give stronger vertical velocity 
and upper tropospheric zonal winds. This is a result of 
stronger cumulusS convection and precipitation whose 
latent heat release provides more See input to the 
upper troposphere. At lower levels, the easteries weaken 
as a consequence of the momentum dissipation at the 
surface. 

When the vertical velocity is intense, a multi-band 
structure often appears. It is characterized by more 
than one concentrated updraft region, stacked closely to 
each other along the ekaudaed slanted frontal surface. 
The stronger updraft is situated equatorward and has a 
larger vertical extent. 

The frontal zone develops equatorward of the dilation 
axis of the large-scale deformation forcing. In the 
midlatitude cases, the frontal zone is closer to this 
axis (about 4° latitude separation at the surface), 
while in the subtropical cases the separation is greater 
(about 5-6°) at the surface. This 1S apparently a result 
of the smaller f and, therefore, a relatively larger 
ageostrophic wind which advects the surface front 


farther equatorward. In addition, the effects of the 


8] 










convection-induced secondary Cc igemhat aon may also 

contribute to the more equatorward position of 

surface front. 

In the ensuing section these numerical results will 
discussed in relationship to the development of the low- 


level westerly jet. 
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- Fig. 3.1 


% Pee " : 
S M60 M90 S60 SSO | 


300.0 320.0 340.0 360.0 380.0 400.0 


Initial potential temperature (solid), 6 (°K), and equivalent 
potential temperature (dash), 6e(°K), soundings at the center 
of the model domain for midlatidude cases - M60 and M90 
and subtropical cases - S60 and S90. (see text). 
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As in Fig. 3.22, except for vertical pressure velocity,-a(10-* cb s7!), 


Fig. 3.27 
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Fig. 3,32 Surface fluxes of momentum (dash; -107* m cb s7?), 


sensible heat (solid; -10-4 °K cb s74) and latent 


heat (dot; -10-4 °K cb s-4). (a) S61 and (b) S91. 
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IV. MECHANISM FOR THE LOW LEVEL JET (LLJ) 


The previous section shows that a low-level westerly 
regime tends to develop in the frontal zone due to Gite 
intensifying meridional temperature gradient in the 
lowest few kilometers. In many of the subtropical cases 
this low-level westerly regime actually develops ina a 
low-level jet whose height of 2-4 km is an good agreement 
with the levels (800-600 mb) of the observed low-level 
jet @Chen, 977; Taig): The jet maximum of 10-15 m/s in 
these subtropical cases is also in reasonable agreement 
with observations. 

In the Introduction, it was pointed out that there are 
several theories for the development of the low-level jet 
nl ae Since effects of terrain (Holton, 1967; Wexler, 
1961) and diurnal cycles (Blackadar, 1957) are not 
included in the model, the related boundary-layer 
theories cannot be treated in this study. Neither can 
the theory of the downward transport of momentum by 
cumulus convection (Masomoto and Ninomiya, 1969; 
Ninomiya, 1971), because the convective adjustment used 
in the present model does not include the 
parameterization of cumulus momentum transport. 

The theory proposed by Uccellini and Johnson (1979) in 
the study of midlatitude severe storms over the U. S. is 


based on a wind-mass adjustment, whereby the secondary 
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circulation forced by an upper level jet streak causes a 
vertical coupling, which leads to the development of a 
LLJ. In this mechanism, the LLJ will develop beneath the 
upper level jet streak such that the two jets will be 
aligned vertically and the axes will be perpendicular to 
each other. Although the model-produced low-level 
westerly maximum often appears (especially in the 
midlatitude cases) as a downward extension of the upper. 
tropospheric westerly jet, the upper level jet is always 
situated much too poleward of the LLJ core for wind-mass 
adjustment mechanism to be effective. 

In the present study, the LLJ in the subtropical 
experiments is much better defined than in the 
midlatitude experiments. Furthermore, in all cases when 


the vertical motion is enhanced due to increased relative 


humidity (therefore also enhanced conditional 
mas tability.) and/or surface fluxes, the low-level 
westerlies become stronger. This is also true with the 


development of the upper level easterlies equatorward of 
the front, which also appears with small magnitudes in 
the midlatitude experiments with higher melative 
humidity . In all subtropical experiments, these upper- 
level easterlies expand and dominate most of the upper 
troposphere equatorward of the front. These results 
suggest that when the vertical motion above frontal zone 


is intense, a “reversed Hadley cell" exists equatorward 
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of the front, with rising Motion Bae newt conw freee 
and sinking motion further equatorward. such an 
overturning would provide an easterly Coriolis torque at 
the upper-level equatorward flow and a westerly Coriolis 
torque at the lower-level poleward flow, which may 
explain the upper-level easterlies and the low-level 
westerly jet. Indirect evidence of the "reversed Hadley 
cell” to the south of ELS may be found in an 
observational study by Ninomiya and Akiyama (1974). They 
found that moist air is located above the LLJ in the 
middle and upper troposphere, but the air south of the 
LLJ is fairly dry throughout the whole troposphere. This 
is indicative of a large-scale subsidence south of the 
PY o Bue. 

More evidence of this "reversed Hadley cell" is seen 
in the temperature distribution of all cases in which 
strong cumulus convection is indicated. For example, the 
temperature in Fig. 3.18 (S90 - the typical subtropical 
case) increases equatorward throughout the troposphere 
for the first day. Afterward, reversed meridional 
temperature gradients begin to appear in the tropical 
region, such that by day 5 the maximum temperature at any 
level in the middle-upper troposphere is near the frontal 
zone. This situation is also apparent in Fig. 3.25 for 
all subtropical cases. The reversed temperature gradient 


to the south of the maximum upward motion has been 
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observed both in the Mei-~Yu front (Chen and Chang, 1980) 
and in the Baiu front (Matsumoto et al., 1970; Ninomiya 
and Akiyama, 1970). Therefore, the equatorside "reserved 
Hadley cell" is actually thermally direct, with warm eu 
rising in the frontal zone and cool air sinking 
equatorward. Hence, this circulation can be maintained 
and as well as continue to provide the upper-level 
easterly and lower-level cuit Pe va cionee 

The next question concerns the identification of the 

main mechanism(s) for this equatorward overturning. 
There appears to be three possible mechanisms: 

1) The secondary circulation normally associated with 
the large-scale flow of a front; 

2) The cumulus convection due to the conditional 
instability, which can increase the latent heat 
release and HOree an enhanced secondary 
circulation; 

3) The third possibility is based on the observation 
that the vertical velocity field of several of the 
more intense cases exhibits a Men ce an a 
structure. Each of the rising branches 1s 
upright, and is stacked in a poleward slanting 
configuration. In a way this structure resembles 
the "slantwise convection" associated with 


symmetric instability in a baroclinic atmosphere 
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(e.g., Emanuel, 1985). This slantwise convection 
may also be a relevant mechanisn. 

To differentiate the last two mechanisms, the term 
“upright convection” will be used to describe the 
basically single-band, vertically-oriented convec tii 
which is normally associated with conditional 
instability, and Cee Wiisic convection" will be used to 
describe the slanted and stacked multi-band structure 
which resembles the convection predicted by the symmetric 
instability theorv. Fig. 4.1 illustrates schematically 
the differences among the three mechanisms. 

The meridional-vertical streamfunction 1s shown in 
Figs. 4.2 and 4.3, for the midlatitude and subtropical 
cases, respectively. For the midlatitude R.H.=60% cases 
(M60 and M61), the atmosphere remains statically stable, 
and Figs. 4.2a,c indicate that the large-scale frontal 
secondary circulation consists of a single cell in the 
Hadley sense. When R.H.=90%, the initial sounding of the 
midlatitude cases (M90 and M91) is conditionally unstable 
in the lower latitudes. At day 4, both cases show a 
small and very weak "reversed Hadley cell" in the middle 
troposphere equatorward of the front. Apparently this 
cell is responsible for the small region of upper level 
easterlies equatorward of the front. However, the modest 


cumulus convection through the release of conditional 
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instability is not intense enough ‘to Pioneer a Sstron zg 
reversed Hadley cell that can develop a strong LLJ. 


All the subtropical cases are initially conditionally 


unstable. A larger and stronger "reversed Hadley cell” 
is present in all cases (Fig. 4.3). This celery is 
Situated poleward and 1S more vertlca Miyerconf ined ™tfor 


the R.H.=60% case than for the R.H.=90% case with the 
middle-level rising motion eein poleward of the 
sumiace front. This is the situation when the LLJ is 
either ill-defined (S60 case, Fig. 3.28a) or rather weak 
(S61, Fig. 28c). When R.H.=90%, the reversed Hadley cell 
is located more equatorward, is more intense and extends 
throughout the verical domain (Figs. 4.3b,d). The LLJ 
is well-defined in this situation (Figs. 3.28b,d), and 
clearly coincides with the eed poleward flow a 
the lower branch of the reversed Hadley cell. The 
cumulus convection is the strongest in these cases and is 
apparently responsible For generating the intense 
equatorside meridional circulation which, in turn, leads 
to the development of the LLJ. 

Ist 1s clear that the Cross Tronta L secondary 
circulation is different from the secondary circulation 


equatorward of the front, and therefore it can not induce 


the LLJ. However, the roles of upright convection and 
slantwise convection are not easy to distinguish. The 
slantwise convection structure 1s present in both 


eel 


midlatitude (when the LLJ is not present) and subtropical 
(when the LLJ is developed) cases. The slanting 1s 
apparently related to the structure of the # frontal 
surface. Although the frontal forcing cannot, by itself, 
directly cause the equatorward secondary circulation, it 
exerts a strong influence on the environmental conditions 
that affect convection. This may involve two mechanisms. 
The first is the large-scale ieee of moist air. This 
leads to either stable precipitation or a "C Loe 
(Conditional Instability of the Second Kind) type 
convection, where the low-level convergence is provided 
by the frontal forcing. The second involves the creation 
of conditions for symmetric instability. These conditions 
occur when the vertical gradient of equivalent potential 
temperature ( @6@e) becomes nesatiae along the Sire One 
pseudo-~angular momentum, which 1s defined by the 


following function (Emanuel, 1983a,b): 
Me=u + fy. 


In a symmetrically unstable baroclinic atmosphere, 
slantwise convection occurs when there is a slantwise 
displacement of an air parcel along a constant M surface. 
Based on a diagnostic semi-geostrophic model, Emanuel 
(1985) further showed that even in the presence of small 
symmetric stability, frontal forcing can result in a 


strong, concentrated, sloping updraft which is_ located 
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ahead of the region of maximum geostrophic fronto- 
genetical forcing. This possibility is supported by 
Sanders and Bosart (1985) who reported observations of a 
concentrated sloping updraft that resembles slant~-wise 
convection in a mesoscale precipitation band. The sloping 
updraft remains even after the environment becomes 
Symmetrically stablized. Apparently the fronto-genetical 
process was force the moist pee nitt “Signed the sloped M- 
surfaces where narrow convective bands can develop. 
Thus, slantwise convection may be expected in areas of 
neutral or small symmetric stability when frontogenetical 
forcing is present. 

The evolution of Be and the M-function for the 
typical midlatitude and subtropical cases (M60 and S90) 
is shown in Figs. 4.4 and 4.5. In the midlatitude case, 
the frontogenetical area is statically stable initially 
and the M-surfaces tilt only slightly (Fig. 4.4a). As 
the frontogenetical process proceeds, the frontal 
structure develops and leads to a stronger tilt of the M- 
surfaces in the highly-baroclinic frontal Zone. This 
results in an area of near-neutral symmetric instability 
after day 1, although the area is limited to the lowest 
troposphere. In the subtropical case, conditional 
instability is present throughout the domain initially 
@ia gummed . 5a). This condition continues until day 4 when 


the frontal structure is developed (Fig. 3.19) and the 
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frontal zone becomes stable in the upright. convection 
sense. As a result of frontogenesis, the M-surfaces tilt 
more towards the cold air in the frontal zone. This 
increased tilt leads to a condition of symmetric 


instability in a relatively deep layer between the 


surface and 6 km in the frontal zone. As the convective 
instability 1s eliminated, a slantwise convection 
structure develops in the vertical motion field. Lat 


appears that this structure 1s produced by the 
frontogenesis forcing of moist updrafts along the sloped 
M-surfaces, similiar to Sanders’7 and Bosart’s (1985" 
observations. 

Distribution of day-5 | Be and the M-function are 
presented in Figs. 4.6 and 4.7 for the midlatitude and 
subtropical experiments, respectively. In the R.H.=60% 
midlatitude cases (M60 and M61, Figs. 4.6a,c), static 
stability remains positive and the area favorable for 
slantwise convection is rather limited and is confined to 
the lowest 2 km in the frontal zone. Thus, the secondary 
Circulation is almost entirely a result of the frontal 
Forernd. 

Distribution of Ge and the M-function are given in 
Figs. 4.6b,d for the R.H.=90% midlatitude cases (M90 and 
M91). In these cases, the equatorward side of the ?eam: 
is conditionally unstable while the poleward side is 


stable initially. Tt ais iunteresting to note that in the 
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frontal zone the convective adjustment apparently has 
stablized the atmosphere in the vertical direction. 
However, the §8e and M lines are nearly parallel in the 
frontal zone, and the area favorable for sloping updraft 
is much greater than that for the R.H.=60% cases. In 
this case, a low level maximum wind of 5 m/s actually 
develops near z=3 km. 

Similar to the a “90s aa@llatitude cases, all Te 
subtropical without-flux cases (Fig. a 7.,.D) have 
conditional unstable profiles equatorward of the front 
and a stable Ge profile poleward. In the with-flux cases 
(Fig. 4.7c,d), the conditional instability also exists in 
the poleward region, a result of the supply of moisture 
and heat at the surface. | The frontal rising motion zones 
oe all cases again —— stable 6 .e profiles, but shows 
signs of even stronger symmetric instability in lower 
troposphere. 

In the S60 case (Fig. 4.7a), there is a larger area of 
near neutral 6e profiles along the M-surfaces in the 
vicinity of the frontal rising motion. In the S6l case 
(Fig. 4.7c), the @e profiles in the lower troposphere 
actually become unstable along the M surfaces in the 
rising motion zone. The same is true for the S90 and S9l 
cases (Figs. 4.7b,d), except the area of symmetric 
instability is even larger than in the R.H.=60% cases. 


Thus, for all cases when the condition of symmetric 
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instability is suggested by the 6e profile along ag 
surfaces, a relatively intense LLJ is developed. While 
in the present model framework, it is difficult to assess 
directly the possible roles of symmetric instability, the 
appearance of the slantwise convection structure seems to 
be an important factor in the occurrence of the 


equatorward secondary circulation. 


Ez 














Fig. 4.1 Schematic diagram to show the differences between 
(a) frontal secondary circulation, (b) upright convection, 
and (c) slantwise convection. 
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V. SUMMARY AND CONCLUSIONS 


In this work a two-dimensional frontal model is used to 
study the structure and behavior of the Mei-Yu front, with 
special emphasis on the dynamics of the low-level jet (LLJ) 
that is frequently observed ahead of the front. The Mei-Yu 
front develops in the area of the east Asian monsoon trough 
during early ier and ion the most important rain- 
producing system over east Asia outside of the typhoon 
season. This study 1s motivated by the unique mixed 
midlatitude-tropical properties of the Mei-Yu front (e.g., 
Chen and Chang, 1980), where heavy convection and strong 
low-level PS Lisonee shear are associated with the nearly 
zonally-oriented front. The occurrence of the LLJ is highly 
correlated with heavy convective rainfall. The mechanism for 
the development of this LLJ is the main focus of this study. 

In the numerical model, the time-varying response to @ 
constant large-scale forcing was treated by ae set of 
perturbation equations. The large-scale forcing was 
prescribed by aestretching deformation field which is 
characteristic of the frontogenetical forcing of the Mei-Yu 
f romE: The quasi-steady state responses to this forcing 
were obtained by integrating the perturbation equations from 
an initial state which is given by the seasonally-averaged 
zonal flow. Two major sets of experiments were conducted. 


The first is based on midlatitude conditions, where the 
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Coriolis parameter at 46° latitude and an initially stable 
sounding were chosen. The second is based on subtropical 
conditions, where the Coriolis parameter at 22.5° latitude 
and an initial sounding that is conditionally unstable was 
chosen. For each set of initial conditions, a number of 
supplemental cases were integrated to study the effects of 
varying initial humidity values and surface fluxes’ of 


momentum, heat and moisture. In the midlatitude cases, the 


quasi-steady state 1s reached within feo UY days Ot 
integration. In the subtropical cases, small-scale features 


due to intense convection cause the motion field to 
fluctuate throughout the integration. A 24-hour averaging 
was applied at day or’ (96-—)20°h)" This is considered the 
quasi-steady state as further integration shows that the 
motion remains quasi-steady. 

For both midlatitude and subtropical experiments, the 
frontal zone develops equatorward of the dilation axis. The 
Separation between the front and the dilation axis is 
greater for the subtropical cases than for the midlatitude 
cases. This larger separation is apparently due to the 
-aarver i “dvd, thererore, “a relatively™large agéostrophic 
secondary circulation around the frontal zone whose lower 
branch advects the surface front farther equatorward. The 
Deweatitudermiront extends deeply into™the upper troposphere 
with a strong poleward tilt, whereas the subtropical front 


is confined to the lower troposphere with less tilt, which 


7 


is in good agreement with observations (e.g., Chen and 
Chang, 1980). Along the sloping front, slantwise updrafts 
develop with a multi-band structure. This updraft is more 
evident in the subtropical cases and in the more moist 
midlatitude cases. 

A westerly jet in the upper troposphere and an easterly 


jet in the lower troposphere develop in both midlatitude and 


subtropical cases, with smaller intensities in the 
subtropical cases. The inclusion of higher humidity or 
surface fluxes leads Be stronger upper westerlies, 


apparently as a result of stronger cumulus convection. The 
lower level easterlies weaken both in the midlatitudes and 
the subtropics with the inclusion of surface fluxes. This 
is apparently due to the effect of momentum dissipation at 
the surface. An interesting feature in the zonal wind field 
in the subtropical cases is the concurrent development of 


upper-level easterlies and low-level westerlies equatorward 


of the front. The low-level westerly maximum at z = 3-4 km 
resembles a LLJ, whose intensity increases when higher 
moisture is included. In the midlatitude cases, only a 


small downward and equatorward penetration of the upper 
westerlies appears in the similar area. 

While several theories have been proposed for the 
development of low-level jets in the Mei-Yu and other 
circulation systems, the choice of the model set up prevents 


consideration of some of these theories as a_ possible 
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explanation for the simulated LLJ. The absence of terrain 
and diurnal radiation cycles in the model precludes’ the 
boundary-layer theories (Wexler, TI6r ; Holton, eo: 
Blackadar, 1957), and the absence of a cumulus momentum 
transport parameterization precludes the mechanism due to 
cumulus momentum transport (Masomoto and Ninomiya, 1969; 


Ninomiya, Power) In addition, the Vertical coupling 


mechanism suggested by Uccellini and Johnson (1979) also 


does not apply in the Mei-Yu simulation, because the upper 
tropospheric westerly jet is located too far from the 
Watitude of the LLJ (Fig. 1.4). 

The appearance of an easterly regime in the upper 
troposphere on the equator side of the front above the LLJ 
suggests that the LLJ may be the result of a secondary 
circulation that resembles a "reversed Hadley” cell. This 
possibility is supported by the temperature distribution, 
Were indicates that the ™reversed™Hadley “cel? is™actually 
a thermally Grrect c1lreulaviron that is capable of 
maintaining the kinetic energy. The lower branch of this 
circulation is directed poleward toward the front, and it 
exerts a westerly Coriolis torque which can develop and 
maintain the LLJ. A meridional-vertical streamfunction 
reveals the structure of this secondary circulation, with a 
strong lower brad®h return flow coinciding with the 


development of a LLJ in the more moist, subtropical cases. 


Pow 


Several possible mechanisms may generate the reversed 
Hadley cell. The possibility that the reversed Hadley cell 
is a secondary circulation associated with the basic frontal 
structure was ruled out because in the drier midlatitude 
cases the frontal secondary circulation is a single cell 
around the front with upward motion on the equator side and 
downward motion on the pole side. No LLJ is developed in 
these cases. In the more noe caiees the subtropieen 
cases, the LLJ developed with the equator-side secondary 
cise u lation: This suggests the importance of cumulus 
convection in developing the reversed Hadley cell and the 
LLJ. 

lim al cases where a LLJ is developed, there is a 
sloping updraft with multi-band rising motion cells. This 
suggests that, besides regular "upright" convection due to 
conditional instability, slantwise convection arising out of 
symmetric instability may also be relevant. The roles of 
these two types of convection were studied by examining the 
evolution of the equivalent potential temperature 
distribution in the vertical and along the pseudo-angular 
momentum (M-function) surfaces. In the subtropical cases 
where the atmosphere is conditionally unstable initially, 
upright convection develops rapidly and the conditional 
instability is quickly released. During the early stages of 
the development, the M-function surface is nearly vertical, 


so that the condition for symmetric instability, (d8e/dz)m < 
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0, is also given by that for conditional instability, (d8e/ 
Pao Aswethe frontogenetical process continues, the M- 
surfaces in the frontal zone begin to slope and the 
condition for symmetric instability become different from 
pee tor Conditional instability. As the conditionally 
unstable soundings adjusts to neutral, the @e sounding 


along M surfaces is either neutral or slightly unstable in 


the middle-lower troposphere below 6 km. At this stage, 


slantwise convection is observed in the frontal zone and is 
indicative of the possible importance of symmetric 
instability processes. 

In the more Wigels Sak midlatitude cases, slantwise 
convection is also observed. In these cases, the frontal 
zone remains statically stable. However, the frontogenetic 
forcing produces slanted M surfaces there, with the @Ge 
profiles on these surfaces reduced to only slightly stable 
based on the symmetric instability criterion. This result is 
in agreement with the recent theory of Emanuel (1985) who 
showed that, under frontogenetic forcing, sloping updrafts 
can occur in a concentrated region on the warm side of the 
frontal zone, and result in an enhanced secondary 
circulation. It also agrees with the observation by Sanders 
and Bosart (1985) of a mesoscale precipitation band in a 
winter-time extratropical cyclone over northeastern United 


States, in which a sloping updraft was observed in an 


14] 


environment in which the symmetric stability was weakly 
positive. 

The main findings of this study are that the LLJ in a 
subtropical Mei-Yu front can develop - a result of the 
Coriolis torque of a secondary circulation on the equator 
side of the front. This secondary circulation is thermally 


direct and is in the reversed Hadley sense. It is driven by 


intense condensation heating due to heavy precipitation. 


The upward motion takes the form of slantwise convection as 
a result of the frontogenesis forcing. The slantwise 
convection begins 6 appear after the atmosphere is 
stablized by the convective adjustment. It is conceivable 
that, without the slantwise convection, which is 
concentrated along the frontal surface, the latent heat 
release would spread out over a broader area. This would 
give rise to a weaker secondary circulation. These 
conclusions are consistent with the observations that low- 
level jets are often associated with intense convection and 
heavy rainfall in the Mei-Yu front, and the  indiiived: 
evidence of a sinking region south of the Baiu front as 
revealed by Matsumoto’s (1972) moisture analysis. Further 
study is needed to clarify the possible roles of the 
symmetric instability mechanism. However, the strong in- 
phase correlation between the rising velocity and heating 
indicates that the main conclusion concerning the role of 


the secondary circulation should not depend on the details 
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of the formulation of condensation heating in the model. In 
mac rc—m—me tronmteal fomcingmis the dominating factor that 
determines the region of maximum rising motion. Given the 
moist subtropical environment for the Mei-Yu system, if a 
specified (instead of predicted) heating function is 
appropriately prescribed in the model, a strong secondary 


Sireubation- equatorward of the front may also be expected. 


Observational studies (e.g., Chen and Tsay, 1978) 


reported that the Mei-Yu convection and heavy rainfall are 
organized in mesoscale convective complexes (MCC), which can 
not be treated in the present two-dimensional model. A 
three-dimensional mesoscale numerical predication model, 
where the MCC’s can be explicitly simulated, should be used 
Mememecxamine the relationship between the LLJ, secondary 
Gmrmeulation and convection. A three-dimensional model is 
also necessary to study the entire life cycle of the Mei-Yu 
front and the LLJ. The development and decaying stages of 
the LLJ and the Mei-Yu front are of different time scales. 
While the Mei-Yu front usually persists for several days to 
one week or more, the appearance of the LLJ 1s more 
Sete reneimenature. Phe time evolution of the LLJ within 
the quasi-steady Mei-Yu front is therefore an interesting 
topic for further studies. 

The model may also be improved to study additional 
effects. The inclusion of a planetary boundary layer 


parameterization that is more sophisticated than the bulk 
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aerodynamics and is coupled with a diurnal radiation cycle 
will better represent the effects of friction, especially 
the interaction with radiation processes. The inclusion of 
terrain may be particularly important, since the Mei-Yu 
front is over the mountainous region of southern China, 
Taiwan and Japan. The topographical effect of the mountains 
and the land-sea complex in this region undoubtedly play an 
important role in -both the ARSENE front and the 
embedded mesoscale convective systems. In addition, Keyser 
and Anthes (1982, 1986) have argued that a stratification- 
dependent boundary layer effect is important for the 
intensity of vertical velocity in a surface front, while Mak 
and Bannon (1984) argued that condensation heating alone 
plays a more critical role. A sophisticated planetary 
boundary layer parameterization may help to clarify this 
issue. Finally, further observational studies are also 
necessary to investigate the Mei-Yu processes both in the 
planetary boundary layer as well as on the entire 


troposphere. 
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Appendix 


A. SPECIFICATION OF MODEL GRIDS 


The horizontal grids are uniforn. The grid interval is 
60 km. The pressure, p, 1s used as vertical coordinates and 
k is the vertical grid index. The pressure of kth level, 
pr, in the model is specified by the exponential 


relationship with height, ax, 


5 
Pr = P, @ (A.1) 


where po is the pressure at the model botton. 


Apply (A.1) to the total height of the model, 2, 


— a 
ee 





where pt is the pressure at model top. 
Using uniform increment of height in vertical, the 


corresponding pressure 18 








where kt is the total number of vertical levels in the 


model. 
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In this study, the total number of vertical levels, kt, 
homo. The corresponding pressure levels used in the model 


are shown in Table A.l 


TABLE A.1 FORECAST LEVELS (cb) 


LEVEE p(cb) Dalkey) p (ebe 
1 95.86 19 44,81 37 20.95 
2 91.9 20 42.96 38 20.08 
3 88.1 a 41.18 39 191926 
4 84.45 22 39.48 40 18.25 
5 80.96 23 37.85 41 17.69 
6 77.61 24 = «36.28 42 16.96 
7 74.4 25 = -34.78 43 16.26 
g 7S 26 33.34 44 1 Smee 
9 68.37 27 31.96 45 14.94 
10 65.54 28 30.64 46 ny 
11 62.83 29 29.37 47 13.73 
12 60.23 30 28.16 48 13.16 
13 37, Se a 26.99 49 12.62 
14 55.35 32 25.88 50 12a 
iV 53.06 33 24.81 51 11.60 
16 50.87 34 23.78 52 fia 
17 48.76 35 720 53 10.66 
18 46.75 36 21a 54 10.22 





B. GRID ARRANGEMENTS 

The model grids and the all field variables are arranged 
as in Fig. B.l. The verta&cal leved, ip, is the mi diesen t 
between levels of p. Vertical velocity is defined at p 
levels and the other variables u, v, 8, q and 7 are defined 
me p levels: The horizontal grids are defined by y and ¥. 
Variables 8, q, wand ® are defined at y where ®Oj,k is 
defined at y and p coordinates ‘and is calculated from the 


continuity equation. 


i A (Fy =p) 
Dik > eee t ~ —k—__k-1__ (V5 ok-1 - Vj-l,k-1) 
¥j ~ ¥j-1 
Yj-lVj-l YG Yj) yaaa Yad 

Pkt 
Pk+1 
Pk 
Px 
Pk-l 
Pk-l 





Piha B:. ) The model grid arrangement. 
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C. FINITE DIFFERENCE FORMULATION 

The equation set (2.1.16) and (271.19) =(2. 1). 23). cane 
integrated by a pure marching process. In this study Wa 
spatially staggered finite-difference scheme 1s used as 
compared to the unstaggered scheme used in Williams (1974). 
The variables 6, tw and q are carried at the same points in 
y, and u and v are carried equidistant between these points 
(see Fig. B.1). Winninghoff (1968), Arakawa and Lamb (1977) — 
and Schoenstadt (1980) have shown that this staggered grid, 
which is known as_ scheme BB, is much better than the 
unstaggered grid (scheme A) for geostrophic adjustment. 
This scheme is especially superior when small-scale heating 
occurs, SUCH face tear 94 Sieeomu: The nonlinear terms are 
arranged to conserve mean squares in the advection terms. 
Leapfrog time differencing is used, except the model is 
restarted every fifth timestep with an Euler backward step 
following Arakawa and Lamb (1977). 

For the finite difference formulation, we discuss only 
the spatial derivative terms in the u and v equations [{i.e., 


(21 Se anid eZ lee ale 


—_—- ——a=_ = Seta ————— os = 


Sy VFR k = Clvjst jk + Vik) (Fist,k + FG,K) 


~ (vp K + vpn k CFG yk + Fp-1,k)I/4(ygja1 - ¥5) 


148 





0 
Lay OF ]jk = 


C(o j+1,k+1 +5 jk41) (Fp okt + Fj,k) 


= eyed k 45k) (Fy kt Fy po1) 1/4 pKa = DK) 


ei 


Oy 


F A A 
SI; eam ke Fae yaa + 5-7) 


2? F ' (Fi+i,k - Fy,k) 
5 j,k = ae. ee . _ 
Oy (Yj+1 > y;) 
0 OFs. _ K+] (Fj ,k+1 - Fy yk) 7 
1,6 ~ 
P (Pk+1 - Pk) 


where F is either u or v. 


149 


For the 8 and q equations, 
0 a x 
Lay VEIjsk = EVik(Gj+1,k + 85k) - Vj-1,k (Gj yk + Gj-1,k)/2(¥j - ¥j-1) 


a n is 
sy #845 ,k = £541 (95 ken + G5 ok) - OF K(GGok + Gj k-1)/2(Pk41 - Pk) 


: 
Sei sk = (G54. jk - Sj-1,k)/(yjet + ¥j-1) 


moths k - 93 ,k) (G5 4k - Gj-1,k) 
oat ay li, tg Yel - Jj = 5p = pe es yj- 1) 
Ky+y (G5 - Gj.k) Ke (Gs, - Gs 4-7) 
0 0G CL adie Kae! jk k 53 ,k j,k-1 
Sp K Spli.k = “ SS SS | 
(Pk+1 - Pk) (Pk - Pk-1) | 


where G is either 6 or q. 
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